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TECHNICAL MEMORANDUM 


ANALYSIS OF GRAPHITE-REINFORCED CEMENTITIOUS COMPOSITES 

1. THE PROBLEM 


l.I Introduction 

In research on composite beams, Biszick^ noted that the transform section method did not 
accurately predict the behavior of the test articles being considered. Biszick documented that a highly 
efficient stmcture; i.e., weight and strain energy, could be developed using a matrix with a very low 
Young’s modulus (E-^ and reinforcement with a very high Young’s modulus (E2)- These test articles 
were cementitious beams reinforced with a graphite mesh. The graphite mesh was used in the same way 
as steel mesh in conventional reinforced concrete; i.e., the reinforcement reacts on the tensile loads, 
while the concrete reacts on the compressive loads. Significant differences should be noted between a 
graphite mesh and a similar steel mesh. First, the graphite mesh has no bending or compression stiffness. 
The mesh drapes, if not supported, and does not resist a compressive force; i.e., cloth. Second, the 
graphite appears to be saturated by the cement (matrix), such that an adhesive bond forms between the 
cement and graphite. The cementitious mix used in the test articles was also unique and exhibited ductile 
properties. Biszick’s research showed conclusively that maintaining a high elastic modulus ratio (77), 
i.e., £^fibe/^matrix’ between the matrix and the reinforcement results in efficient composite stmctures. It 
follows then, to fully utilize this discovery, the analytical methods must be understood and developed. 

This Technical Memorandum (TM) investigates the analytical methods for the analysis of 
graphite-reinforced cementitious composites where high ty’s are employed. A sample problem will 
demonstrate the transform section method of analysis. A finite element model will then be developed 
and the results compared to the transform section method of analysis. The errors associated with the 
transform section method and high rj’s will be demonstrated. Alternative methods of analysis will then 
be investigated. The rule of mixtures^^ will be applied to a graphite -reinforced cementitious composite 
to calculate effective material properties. These properties will be introduced into the transform section 
method and the results compared to finite element methods and beam testing. Laminated composite plate 
theory will then be discussed and applied to the graphite -reinforced cementitious composite. It will be 
shown that laminated plate theory for composites can successfully be applied to this type of composite. 
Extensive experimentation is used to develop the material properties used in the laminated plate theory. 
Graphite-reinforced cementitious composite beams are analyzed using classical composite methods and 
compared to finite element methods and experimental results. 



1.2 Problem Statement 


Section 1.2 illustrates a problem with the analysis of a composite beam that has a high elastic 
modulus ratio. First, the expression for the moment of inertia (/) is derived in terms of rj. The moment 
of inertia is then used to calculate the deflection of a beam subjected to pure bending. The deflection of 
the beam is calculated as rj ranges from 1 to 1 ,000. The analysis is then repeated using a finite element 
model. The conclusion is made that the transform section method of analysis is not accurate for rj > 20. 

1.3 Composite Beam 

Figure 1 illustrates a beam composed of two different materials. Since the beam is subjected to 
boundary conditions and loads to produce pure bending between the supports (fig. 2), the transform 
section property method may be used to determine the deflection in the beam for each of the load cases 
described in table 1 . These results can be compared to those obtained from a finite element model.^ 
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Figure 1. Beam cross section. 
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Table 1 . Loads and deflections for transform section method. 


Case 



n 

Flange 


Deflection 

1 

3.00E+07 

3.00E+07 

1 

0.347 

0.0289 

9.22722E-03 

2 

1.50E+07 

3.00E+07 

2 

0.694 

0.0542 

9.8401 OE-03 

3 

7.50E+06 

3.00E+07 

4 

1.388 

0.1048 

1.01 781 E-02 

4 

3.00E+06 

3.00E+07 

10 

3.470 

0.2566 

1.03923E-02 

5 

2.14E+06 

3.00E+07 

14 

4.858 

0.3578 

1 .04341 E-02 

6 

1.50E+06 

3.00E+07 

20 

6.940 

0.5097 

1.04657E-02 

7 

1.00E+06 

3.00E+07 

30 

10.410 

0.7627 

1.04904E-02 

8 

7.50E+05 

3.00E+07 

40 

13.880 

1.0157 

1.05028E-02 

9 

3.75E+05 

3.00E+07 

80 

27.760 

2.0278 

1.05215E-02 

10 

3.00E+05 

3.00E+07 

100 

34.700 

2.5338 

1.05252E-02 

11 

2.50E+05 

3.00E+07 

120 

41.640 

3.0399 

1.05277E-02 

12 

1.50E+05 

3.00E+07 

200 

69.400 

5.0640 

1.05327E-02 

13 

1.36E+05 

3.00E+07 

220 

76.340 

5.5701 

1.05334E-02 

14 

1.00E+05 

3.00E+07 

300 

104.100 

7.5942 

1.05352E-02 

15 

7.50E+04 

3.00E+07 

400 

138.800 

10.1244 

1.05364E-02 

16 

6.00E+04 

3.00E+07 

500 

173.500 

12.6547 

1.05372E-02 

17 

5.00E+04 

3.00E+07 

600 

208.200 

15.1849 

1.05377E-02 

18 

3.00E+04 

3.00E+07 

1,000 

347.000 

25.3057 

1.05387E-02 


1.4 Analysis Using the Transform Section Method 

The equation for the deflection ^ beam (fig. 3) in pure bending can be derived using 

the elastic curve equations and the boundary conditions: 


^’max 


Pax 

2EiI 


(L-x) ■ 


Force 




Figure 3. Free-body diagram. 




The effective value of the moment of inertia can be calculated using the transform section 
method;^ and A3 are made from a material with properties of £’2- A2 is made from a material with 
properties of E^. The original cross section of the beam (fig. 4 ) is modified per the transform section 
method. The moment of inertia is then calculated from the new geometry and summarized below: 

• Material properties of the composite beam: £j = 360 x 10^ Ibf/in.^, £2 = 30 x 10^ Ibf/in.^ 

• Transform section ratio: 77 = £2/£j, rj = 83.33 

• Width of the composite beam: [i = 0.347 in. 

• Moment of inertia calculation: 



= 3.645x10 ’^ 1773 = 1.0416 xlO’^yS = 3.645 xl0~^ 77^8 

II = 1 .045 iu.^ 1 2 = 0.004 iu."^ = 1 .045 in.^ 


/ = /l + /2 + /3 / =2.112^."^ 



Figure 4 . Transform section geometry. 


It should be noted that the middle section of the beam (A2) is independent of 77. As 77 increases, 
section A2 becomes insignificant. The expression for I can be written in terms of rj as follows: 


4 



I = 3.645 X 10-277^3 + 1.0416 X IQ-^yS + 3.645 x 


/ = 7.29x10-27773 + 1.0416 x10-2;8 . 


( 2 ) 


Let P = 0.347 in., then 


/ = 2.53x10-277 + 3 . 6 x 10-3 in .4 

Now this expression can be substituted into the elastic curve equation: 


Pax 


'max 


2^1 ( 0 . 0025377 + 0.0036) 






( 3 ) 


( 4 ) 


Since 



Pax 


2^1 (0.00253(E2/Ei) + 0.0036) 

v2. 


( 5 ) 


so that 


^max 


Pax 

0.0506^2 +0-0072^i 


v2. ■ 


( 6 ) 


Equation (6) is used in a Microsoft® Excel spreadsheet to calculate the deflection for each case. 
As shown in table 1, £2 remains constant while decreases. As the Young’s modulus ratio (77) 
increases from 1 to 1,000, the deflection value increases to 0.015 in. and then becomes essentially 
constant. 


1.5 Analysis Using the Finite Element Method 

A finite element model of the beam was generated using the MSC/Nastran finite element code.^ 
The model is composed of quad4 plate elements with loads and constraints (fig. 5). The model returned 
a maximum deflection of 9.19 x 10^3 for 77 = 1. The results of the finite element model can be compared 
to the results from the elastic curve equation for the same case: 


Error% = 


Y —Y 
model equation 


model 


xl00% = 


0.00919-0.00922 


0.00919 


X 100% = 0.39% 


( 7 ) 


5 



9.19-03 



Default Deformation: 
Max 9.1 9-03 @Nd 103 


Figure 5. Finite element model. 


The analysis of the first case indicates that the transform section method and the finite element 
model are in good agreement. The finite element model was run for each of the 18 cases. The results 
were tabulated in the Excel spreadsheet (table 2). 

Table 2. Results of the finite element analysis. 


Case 

mm 

E2 

n 


Deflection 

Model 

m 

1 

3.00E+07 

3.00E+07 

1 

0.0289 

9.22722E-03 

9.19E-03 


2 

1.50E+07 

3.00E+07 

2 

0.0542 

9.8401 OE-03 

9.78E-03 

mSm 

3 

7.50E+06 

3.00E+07 

4 

0.1048 

1.01 781 E-02 

1.01 E-02 

-1.10 

4 

3.00E+06 

3.00E+07 

10 

0.2566 

1.03923E-02 

1.02E-02 

-2.01 

5 

2.14E+06 

3.00E+07 

14 

0.3578 

1 .04341 E-02 

1.02E-02 

-2.30 

6 

1.50E+06 

3.00E+07 

20 

0.5097 

1.04657E-02 

1.02E-02 

-2.24 

7 

1.00E+06 

3.00E+07 

30 

0.7627 

1.04904E-02 

1.04E-02 

-1.29 

8 

7.50E+05 

3.00E+07 

40 

1.0157 

1.05028E-02 

1.06E-02 

0.51 

9 

3.75E+05 

3.00E+07 

80 

2.0278 

1.05215E-02 

1.19E-02 

11.88 

10 

3.00E+05 

3.00E+07 

100 

2.5338 

1.05252E-02 

1.29E-02 

18.16 

11 

2.50E+05 

3.00E+07 

120 

3.0399 

1.05277E-02 

1.39E-02 

24.15 

12 

1.50E+05 

3.00E+07 

200 

5.0640 

1.05327E-02 

1.84E-02 

42.76 

13 

1.36E+05 

3.00E+07 

220 

5.5701 

1 .05334E-02 

1.96E-02 

46.26 

14 

1.00E+05 

3.00E+07 

300 

7.5942 

1.05352E-02 

2.42E-02 

56.54 

15 

7.50E+04 

3.00E+07 

400 

10.1244 

1 .05364E-02 

2.96E-02 

64.44 

16 

6.00E+04 

3.00E+07 

500 

12.6547 

1 .05372E-02 

3.44E-02 

69.35 

17 

5.00E+04 

3.00E+07 

600 

15.1849 

1 .05377E-02 

3.85E-02 

72.59 

18 

3.00E+04 

3.00E+07 

1,000 

25.3057 

1.05387E-02 

4.86E-02 

78.31 
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The deflection values from the two methods of analysis are plotted as rj increases from 1 to 80. 
Figure 6 illustrates that the transform sections method of analysis results in a constant deflection as rj 
approaches 20. The finite element model and the transform section method are in reasonable agreement 
until rj reaches 40. Beyond 40 the significant differences in the analysis methods are obvious. 


e- Transform Section — B— Finite Element Model 



Figure 6. Comparison of transform section method to finite element methods. 


This sample problem illustrates the current state of the art for the analysis of graphite-reinforced 
cementitious materials. This analysis leads to the conclusion that the transform section method of analy- 
sis is not adequate for the analysis of graphite -reinforced cementitious composites. The literature search 
included in section 2 was undertaken to determine an alternate approach. 
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2. LITERATURE REVIEW 


Section 2 provides the background necessary to develop an understanding of the process required 
to analyze graphite -reinforced cementitious composites. The basic theories used in the analysis are 
discussed to provide context for subsequent derivation. The material presented here is intended to give 
the reader a sense of direction for this TM. 


2.1 Introduction 

A fiber-reinforced material consists of two components: fiber and cement. The cement fills the 
space around the fibers and keeps them in position. The hardened cement adheres to the fibers, allowing 
a transfer of load. The type of fiber reinforcement may consist of a variety of materials: steel bars, glass 
fabric, ceramic whiskers, or carbon strands, to name a few. There are many types of cement as well, 
including epoxy, polyesters, aluminum alloys, and cementitious materials. In general, the cement mate- 
rial is isotropic while the fibers are longitudinal components that can react to forces along their axes. 
The fibers are typically much stronger and stiffer than the cement. The material properties of the fiber- 
reinforced composite are typically orthotropic in nature. The elastic properties of the composite depend 
on the properties of its constitutive materials, the mix ratio of the components, and the orientation of the 
fibers. 


The literature associated with composite analysis and development covers a phenomenal range 
of applications, from the International Space Station to the common kitchen sink. However, the one 
common and elusive goal of researchers is to develop a method to predict the elastic properties of a 
composite from the properties and geometry of the constitutive materials. The purpose of this literature 
search is to determine which methods will most likely work for the graphite cementious materials being 
considered in this TM. 


2.2 Structural Engineering Approach 

The structural engineering approach is illustrated in the transform section method of analysis 
(sec. 1.3). The most common form of the equation is explained in detail by Beer and Johnston.^ This 
approach is often associated with the analysis of steel bars embedded in a concrete matrix. Allen and 
Haisler^ modify the approach in a concept referred to as the modulus weighted section properties. In this 
approach, the beam is divided into discrete homogeneous parts (fig. 7). The modulus weighted section 
properties are then defined in terms of the Young’s modulus ratios as follows: 


n p 




( 8 ) 


8 



( 9 ) 



Figure 7 . Composite cross section for weight modulus properties. 


This method of determining the section properties for a composite beam is considerably more 
complex than methods described by Beer and Johnston.^ However, the elastic modulus ratio (rj = E2/E1) 
governs this equation, and as the value of rj increases, the accuracy of the method decreases. 

Balaguru^ developed another interesting approach to the analysis of composite beams. 
Ferrocement is a composite material, consisting primarily of several layers of wire or fiberglass mesh 
embedded in a cementitious matrix. Balaguru^ developed an analytical model of the composite beams 
from the stress-strain curves of the constitutive materials. The model uses section properties calculated 
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from the transform section method. However, Balaguru takes into account the nonlinear properties of the 
materials. As the weaker of the materials becomes nonlinear. Young’s modulus is recalculated. The new 
value is used to determine the section properties for the nonlinear load. Balaguru ’s methods seemed 
promising for the analysis of the composite considered here; particularly, since the method had been 
used for the analysis of glass fibers with a cementitious matrix. However, on further investigation,^ it 
was found that the analytical method suffered from the same limitations as those encountered with the 
transform section method. Nonetheless, Balaguru ’s methods were successfully used with a modified 
transform section method (sec. 4.13) to determine the deflections of a composite beam considered here. 

The common theory used for the structural engineering approach is the transform section 
method. This approach is the primary cause of “the problem” discussed in section 1. Further investiga- 
tion into the structural engineering approach would be unavailing. 

2.3 Mechanics of Materials Approach to Composites 

The mechanics of materials approach to composites is fundamentally different from the 
structure’s approach discussed previously. This approach is based on the mix ratio (volume fraction) 
between the fiber and cementitious materials. The mechanics of materials approach is often referred to 
as the rule of mixtures for composites. The goal of this theory is to determine the “effective” elastic 
material properties: Young’s modulus {E^ and £’ 2 ), shear modulus {G-^ 2 ), and Poisson’s ratio (v^ 2 ) 
the material properties of the fiber and cementitious material. Numerous examples for the mechanics of 
materials approach can be found in the literature,^’^*^ and the method is explained in detail in section 4. 
The literature is unanimous^’^*^® in declaring this method accurate for determining the effective material 
property E^, provided the fibers are continuous and aligned in the direction of the load. The remaining 
material properties, E 2 , G 12 , and v^ 2 » be accurately predicted by the rule of mixtures. Although 

Jones^° provides detailed and convincing derivations for these material properties, he points out several 
assumptions in the theory that would lead to errors in the calculation of the values. In addition, Jones^*^ 
strongly points to the need for experimental data to verify the results. Balaguru^ provides the equations 
to determine the value of however, he chose not to address the remaining properties. The literature 
contains numerous variations on the rule of mixtures theory. These variations attempt to correct the 
deficiencies in the theory. 

Halpin-Tsai,^’ll for example, made significant contributions to the rule of mixture methods by 
developing the Halpin-Tsai equations: 


p ~ pr Y + F V 

^ fiber fiber ''' '^matrix' matrix 


^12 — ^fiber^fiber '^matrix^matrix 


£2 _ (l + g^^fiber) 

^ matrix (1 ^^fiber) 


, where 77 = 




"fiber 


V ^matrix J 


" fiber 




matrix J 


( 12 ) 

(13) 

(14) 
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(15) 


^12 _ (l + C^^fiber) 

^matrix “ ^^fiber ) 


where rj = 


^fiber 
V '^matrix 



yv 


^fiber 

r< 

'^matrix 


+ Q • 
y 


The measure of reinforcement (g) depends on the boundary conditions within the composite and 
is an empirically derived value for the specific composite. Jones provided sufficient evidence that these 
equations are accurate, provided g is available for the composite being considered. However, no data 
were available for graphite-reinforced cementitious composites. 

Another promising variation on the mle of mixtures method is provided by Krenchel.^ Krenchel 
used an equilibrium approach that considered reinforcements in both the v and y directions. The methods 
provided by Krenchel were applied to the composite considered here. The ply angles y/i and y/2 are 
measured from the direction of the load. The P term represents the volume fraction for the fiber. The k 
term is given by k = \- P ! (\- Vcomposite ’ must be estimated since is not known. The 

following equations are then solved simultaneously to determine ^^composite- 

• Stress in direction of force -x: 


^composite^x ^matrix^(^x '^matrix^y )+ ^fiber/^cos(v^i)n e^cos()/ri)^ +e^ sin(v^i)^ 


' 2 /f 2 "" 

£fiber^cos(vA2) |^e^cos(vA2) +eySm(Y2y 


(16) 


Stress in transverse direction -y; 


0 ^matrix^('^matrix^x ^fiber^^i^(v^l) +£ySin(l/r]^) 


-I- 


2 / 2 ^ 
^fiber/^cos(v^2) £xCos(vr2) -i-e. sin (1/^2)' 


(17) 


Poisson’s ratio: 


composite g 


(18) 
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The Krenchel approach provides significant promise because it considers the loads and forces 
in both directions. His equations were solved for the example considered in section 4.6 and returned 
correct values for and £'2 (app. A). However, the shear modulus (G 12 ) returned the same value as the 
previous methods. Additionally, Krenchel had no verification data for the shear modulus calculations. 

In summary, it can be said that the mechanics of materials approach will provide a set of effec- 
tive material properties for a fiber-reinforced composite. Balaguru^ and KrencheP applied the method 
for fiber-reinforced cementious materials. Krenchel’s methods can be manipulated to provide Ei and £2 
for a composite. However, the literature indicates that estimates of the shear modulus and Poisson’s ratio 
were questionable. Researchers were unanimous in their endorsement of establishing an experimental 
database for any composite being characterized for the first time. 

2.4 Experimental Methods 

The literature search for experimental methods concentrated on the standard test methods pub- 
lished by the American Society for Testing and Materials (ASTM). Two ASTM standards described the 
possible geometry of a tensile test specimen. ASTM Standard D 3039^ recommended a tensile test 
specimen with bonded tabs while ASTM Standard established the dogbone geometry. The purpose 
of both configurations is to provide a strong grip area, causing the specimen to strain and break in the 
center. Bonding grip tabs to the cementitious specimens were found to be difficult and impractical. The 
grips used on the tensile test machine limited the thickness of the tabs. One specimen with tab ends was 
successfully tested. However, the machined dogbone specimens (sec. 3.3) were much easier to produce 
and worked well. The ASTM recommendation to use extensometers to measure the strain was a signifi- 
cant impact to the tensile test. The extensometers were found to be unreliable because the knife blades 
would not grip the test article properly. Strain gauges provided a much more reliable method for strain 
measurement (sec. 3.5). Strain gauge installation was achieved by following the recommendations 
of Measurement Group Bulletin 309D.^^ In addition. Measurement Group Tech Tip TT-61 1^^ was 
reviewed for applicable information. Tech Tip TT-61 1 appears more suitable for large concrete 
structures. 

Several experimental procedures for determining the shear modulus were offered by Jones. 

The most compelling was producing a specimen with the fibers running at a 45° angle, and testing the 
specimen as a tensile test. Jones^® provides the necessary equations to determine the shear modulus. 
However, ASTM standard D 5379-93^^ recommended an losipescu shear test for determining the shear 
modulus. On further investigations, it was found that the ASTM standard also had problems. Conant and 
Odom^^ identified numerous problems with the ASTM/Wyoming losipescu test fixture. The most 
significant was the tendency for the specimen to twist during testing. They made improvements to the 
design and termed it the Idaho shear test fixture. Ifju^^ describes a serious inadequacy in the instrumen- 
tation recommended in the ASTM standard.!^ Ifju^^ demonstrated that by using a strain gauge rosette, 
which covered the entire shear area, the accuracy of the test could be ensured. In addition, Ifju showed 
that by placing these gauges on each side of the specimen, the effects of twisting could be removed. 

It was concluded from the literature search that the ASTM standard tensile^’^^ and shear test^^ 
for a composite material were the most desirable methods for determining the material properties 
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for a composite material. Interestingly, no reference was found to these tests being applied to graphite- 
based cementitious composites. 


2.5 Laminated Plate Theory 

Halpin, Jones, and Tsai^^^^ each considered the determination of the composite material proper- 
ties as an intermediate step in the analysis of a composite. The next logical step is to use these effective 
material properties to develop a laminated composite plate. Gurdal^^ pointed out that a significant 
advantage of a composite material is the ability to optimize the geometry of the fibers for a given load- 
ing condition. Krenchel’s^ equations provided a crude means to rotate fibers in a cementitious matrix. 
However, his work was limited in that it assumes each layer will have the same geometry. The literature 
search failed to find an example of the laminated composite plate theory being applied to a fiber- 
reinforced cementitious composite. 

The laminated plate theories provided by Jones et al.^’ih20 contained similar derivations 
and arrived at the same equations for determining the effective material properties of a laminated plate. 
However, these authors each chose different methods to solve this complex system of equations. As 
a result, there is a significant difference in the final effective properties determined by the authors. 
Because of these conflicting results, the laminated plate theory is investigated. The equations are 
arranged in such a way as to allow an exact solution to be determined. In addition, Jones et al. developed 
separate solutions, depending on whether the laminated plate was symmetric or nonsymmetric in nature. 
It is not necessary to make that distinction for the equations developed here. As shown by Nettles,^® the 
laminated composite plate theory can be used to determine the deflections, stresses, and strains for the 
multilayered composite plates. However, for complex loads or geometry, these methods are difficult 
to apply. The finite element methods offer an opportunity to analyze complex layups and geometry. 

2.6 Finite Element Modeling 

Matthews^^ illustrated that finite element analysis is applicable to composite materials, provided 
the user adequately represents the properties and layup of the plates. MSC/Nastran^^ was chosen as the 
finite element code because of the layered composite element capabilities. The example problem pro- 
vided by Nettles^® was analyzed using the finite element model. MSC/Nastran defined the required 
material input and the equations from laminated plate theory were used to generate the input data 
for the finite element model. Although the literature contains references to finite element modeling 
of fiber-reinforced cementitious beams, there was no reference for a multilayered plate. 
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3. MATERIAL PROPERTIES TESTING 


The process of preparing and testing the tensile and shear specimens is described in section 3, 
along with details of mixing and placing the cementitious materials. The specimen geometry is briefly 
discussed and presented. The process of producing a graphite-reinforced composite plate is then 
described, followed by a unique method of machining the plates into tensile and shear specimens. The 
use of several types of extensometers and strain gauges is demonstrated for the tensile test specimens. 
A specialized losipescu strain gauge is described for the shear test. The tensile test procedure and the 
results from six test specimens are provided. The shear test procedure and the results for two shear test 
specimens are also included. 

3.1 Mixing and Placing Cementitions Composites 

Mixing and placing the cementitious matrix material is a very important aspect of specimen 
development. The cementitious mix used in this investigation was developed for use in the American 
Society of Civil Engineers’ concrete canoe competition. The mix development is not a part of this 
investigation. Table 3 defines the mix used in all specimens. 


Table 3. Ingredients of the cementitious mixture. 


Cementitious 

Mixture 

Grams 

K-25 Microbubbles 

218 

Portland cement 

558 

Latex 

310 

Fortifier 

150 

Denatured water 

-620 


The constituents of the mix were each weighed prior to the mixing process (fig. 8). The K-25 
microbubbles were then added to the cement in small amounts and stiiTed by hand. Care was taken not 
to break the microbubbles. The mixing process was judged complete by its color. The latex, fortifier, 
and half of the water were then poured together repeatedly to ensure a thorough mixing of the liquid 
ingredients — a very important step. It was found that the consistency of the cementitious matrix was 
smooth and creamy when the latex and fortifier were mixed with water before they were added to the 
dry ingredients. If the latex and fortifier are added separately, the consistency of the mixture is lumpy 
and coarse. 
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Figure 8. Tools used for mixing and weighing ingredients. 


The liquid ingredients are then mixed with the cement and microbubbles. The remaining dena- 
tured water (~50 percent of the original weight) is added until the mixture reaches the proper consis- 
tency. The cementitious mixture should be smooth and creamy with the consistency of warm butter 
(fig. 9). A mixture with too much water will leak out of the mold, leaving voids in the specimen. A 
mixture that is too dry will not “flow” around the graphite mesh, weakening the bond between the 
matrix and the reinforcement. In addition, a dry mixture will tend to displace the graphite as the mixture 
is pushed through the mesh. 



Figure 9. Cementitious mixture with proper consistency. 
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Numerous attempts were made to remove air from the final mix before it was poured. Air in the 
mix was apparent by the high porosity seen in some of the specimens during finishing and sanding (fig. 
10). The mix was placed in a vacuum bottle with a small vacuum applied (fig. 8). This process should be 
considered marginally successful. Too much vacuum pressure caused the mix to splatter and spill in the 
vacuum bottle. Reducing the vacuum to levels that ehminated the splatter had little effect on the amount 
of air in the mix. The mix was placed on vibrating and shaking tables, causing air bubbles to rise to the 
top of the mix. However, it tended to separate the mix (liquids came to the top) such that some remixing 
was necessary before placement. 



Figure 10. Specimen with high porosity. 

Experience indicated that the best opportunity for reducing the specimen porosity was during 
the placement process. The cementitious mix was added to the mold in small amounts. The thin layers 
of cementitious mix were slowly “worked” into the molds and graphite mesh. Working or floating the 
cementitious mix caused the air bubbles to rise to the top, and porosity was eventually reduced 
to acceptable levels. 


3.2 Tensile and Shear Specimens 

Tensile testing is used to determine the composite material properties E^, E 2 , and v^ 2 - These 
properties are needed to complete the compliance matrix (sec. 5.3) for an orthotropic material: 


S = 




0 


0 




(19) 


The first step in the process is to determine the size and shape of a tensile test specimen for an 
orthotropic composite material. The ASTM provides standards for tensile testing for a wide variety of 
materials. A search of the ASTM database^^ provided a number of testing standards for materials similar 
to the considered composite. The mutual goal of these test methods is to develop a test specimen that 
provides tensile strain data, modulus of elasticity, and Poisson’s ratio. ASTM Standard E8-99, “Standard 
Test Methods for Tension Testing of Metallic Materials,”^^ establishes the dogbone configuration for the 
tensile test specimen, as shown in figure 1 1 . 
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Figure 1 1 . Dogbone configuration. 


This test method employs small specimens that have a reduced cross section at the center of then- 
length to avoid failure in the grip area. The transition from the full width of specimen to a reduced 
section at the center is gradual to minimize stress concentrations. ASTM Standard D 638-99, “Standard 
Test Method for Tensile Properties of Plastics,”^ demonstrates the use of the dogbone specimen for 
reinforced composites, including orthotropic laminates. A unique characteristic of this composite is the 
relatively wide spacing of the graphite tows. The mesh contains eight yams per inch in both axial and 
transverse directions (8 in. by 8 in. weave). Reducing the width of the gauge section limits the number 
of graphite yams in the tensile test. A decision to hold the gauge width to 1 in. made the specimen 
geometry dimensions similar to the type III specimen of ASTM Standard D 638. In each specimen 
geometry, the configuration is determined by the desire to produce failure in the gauge section without 
inducing stress concentrations in the gauge section. The specimen length is substantially longer than the 
width to minimize the bending stresses due to minor grip eccentricities. The final specimen configura- 
tion is analogous to the specimens described in the ASTM standards,^’^^ with the exception that the 
gauge width was 1 in. 

The geometry of the shear specimen is established in ASTM D 5379-93, “Standard Test Method 
for Shear Properties of Composite Materials by the V-Notched Beam Methods.”^^ This test is often 
referred to as the losipescu shear test. Figure 12 illustrates the geometry of the shear test specimen. 



Figure 12. Shear test specimen. 
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3.3 Tensile Specimen Production 


Specimens were made in a plexiglass mold designed to place the graphite mesh at the center of 
the laminated composite. As shown in figure 13, the graphite was placed in tension to prevent the weight 
of the cementitious mixture from pushing the graphite toward the bottom of the mold. 



Figure 13. Plexiglas mold with graphite tensioning device. 


The cementitious mix was floated to the top level of the bottom spacer and a sheet of graphite 
mesh was placed over the bottom spacer. The graphite tensioning device was then installed. The top 
spacer was added and cementitious mix was floated to the top of the spacer. The tensioned graphite mesh 
was trapped between two layers of the cementitious mix. The composite was allowed to cure at room 
temperature for at least 21 days before the mold and tensioning device was removed. Even though mold 
release agents were used, it was necessary to twist and flex the plexiglass to remove the specimens 
without breaking them. After the specimens were removed from the molds, the graphite was trimmed. 
The plexiglass side of the specimen had a smooth finish, while the top of each specimen had a coarser 
texture. A bandsaw was then used to cut the “cookie-sheet” specimen into 2-in. by 14-in. strips for 
machining. 
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An undesirable aspect of the dogbone geometry is that it requires machining of the composite 
laminate. Laminated composites can be difficult to machine;^ delaminations, notches, and uneven 
surfaces can result from inappropriate machining methods. In an effort to fix edges damaged by a 
milling machine, it was found that a bandsaw could be used to cut the cementitious and graphite com- 
posite. The cuts were smooth and notch free. A dogbone specimen was then cut out using the bandsaw. 
Taking the machining process further, it was found that a router would cut a smooth radius on the 
specimen (fig. 14). 



Figure 14. Router used to machine dogbone specimens. 
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The router left the edges smooth and notch free. On further investigation, it was learned that 
a template, shown in figure 15, used with a table-mounted router, provided an excellent method 
of machining the dogbone specimens. A specimen configured with this template is shown in figure 16. 
The specimen has a 5-in. gauge length, a 1-in. gauge section width, and a 14-in. length. 



Figure 15. Template used for machining specimens. 



Figure 16. Machined dogbone specimen. 
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As previously mentioned, the top side of the specimen had a coarse texture. The textured surface 
made it difficult to attach instrumentation. A flat backing plate equipped with fine grit sandpaper, shown 
in figure 17, was used to finish the top surface of the specimen. This process was also used to ensure that 
the graphite mesh was in the geometric center of the specimen. Occasionally, finishing the top surface of 
the specimen resulted in the top layer being thinner than the bottom. In such a case, the bottom layer was 
sanded to the same thickness as the top. Using shimming material as a guide ensured that the top and 
bottom layers were the same thickness (fig. 18). 



Figure 17. Flat plate used to polish surface. 



Figure 18. Edge view of finished specimen. 
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3.4 Shear Specimen Production 


The geometry of the shear specimen was difficult to machine because of the small radius at the 
base of the notch. Designing a router blade with a cutting diameter of 0.05 in. solved the problem. A 
router template was then machined to the dimensions specified in the ASTM standard^^ (fig. 19). Shear 
test specimens could then be machined using a router and table (fig. 20) in the same manner as the 
tension test specimens. 



Figure 19. Shear specimen router template. 



Figure 20. Shear specimen router table. 
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The router machines the shear specimen, leaving the edges smooth and the dimensions within 
specification (fig. 21). The top surface of the specimen requires sanding for a smooth finish and allows 
adjustment of the thickness if required. 



Figure 21. Shear test specimen. 


3.5 Tensile Specimen Instrnmentation 

The instrumentation for the tensile test specimen must simultaneously measure the longitudinal 
and transverse strain. The ASTM standards recommend using extensometers,^ noting that strain gauges 
provide a viable alternative. However, proper techniques for mounting strain gauges are crucial to 
obtaining accurate data.^^ 

Instrumentation of the specimens became a significant challenge to successful testing. Three 
types of extensometers were available during testing. The 2-in. gauge length was used for axial strain 
measurement (fig. 22). 



Figure 22. Extensometer with a 2-in. gauge length. 
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The extensometer was attached to the edge of the test specimen. The data were compared to the 
axial strain gauge data. Plots of the strain gauge and extensometer data for the axial direction were in 
good agreement (fig. 23), indicating that the strain gauge application methods were acceptable. The 1-in. 
gauge length was intended for transverse strain measurement (fig. 24). This extensometer was attached 
to the face of the test specimen. However, the thin specimen did not have adequate stiffness to support 
the weight of the extensometer. The instrument appeared to deflect the specimen, causing the knife 
edges of the extensometer to lose contact with the specimen surface. As a result, the data from the 
transverse extensometer did not correspond to the data from the transverse strain gauge (fig. 25). 


■A — Strain Gauge — Extensometer 



Figure 23. Axial strain gauge and extensometer. 
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Figure 24. Extensometer with a 1-in. gauge length. 



Figure 25. Extensometer mounted on test specimen. 


The biaxial extensometer, as shown in figure 26, was also considered. However, the instrument 
weight deformed the specimen unacceptably. Efforts were made to support the weight of the instrument. 
Unfortunately, the instmment relies on its weight to maintain contact between the knife edges and the 
specimen. It was too difficult to reduce the deformation in the specimen and maintain adequate weight 
on the knife edges of the instrument. The biaxial extensometer was abandoned in place of the previously 
mentioned extensometers and strain gauges. Strain gauges were attached to both sides of the specimen 
(fig. 27). The strain gauges were 90° Tee rosettes, model CEA-06-125UT-120, manufactured by Micro- 
Measurements. The gauges were oriented in the axial and transverse directions. 
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Figure 27. Specimen with strain gauge installed. 


All the strain measurements were made at room temperature (=75 °F). The strain gauge 
installation was completed by the strain gauge manufacturer’s recommendation (fig. 28).^^ 



Figure 28. Strain gauge installation. 
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The specimen was smoothed using fine-grit sandpaper and was free from irregularities and 
contaminates. An acidic cleaner. Prep Conditioner A, was applied to the gauge region. M-Prep Neutral- 
izer 5A followed the acid cleaner. Micro-Measurements M-Bond 200 and a catalyst were used to bond 
the gauges to the specimen. The four strain gauges provided axial and transverse strain data for each 
side of the specimen. This allowed the data to be corrected for any bending. 

The strain gauge installation method worked well providing no porosity or surface “grit” was 
present. As mentioned earlier, the order in which the ingredients are mixed and the manner in which 
the molds are poured significantly affected specimen porosity and texture. Specimens with high porosity 
or a texture that could not be sanded smooth were discarded. 

3.6 Shear Specimen Instrumentation 

ASTM D 537993 defines the instrumentation for the shear test article.^^ This standard specifies 
a strain gauge rosette with two 45° strain gauges be placed in the center of the specimen. The gauges 
are oriented at ±45° to the loading axis (fig. 29). 



Figure 29. Strain gauge location for shear specimen. 


The recommended strain gauge configuration requires the use of correction factors to account for 
the gradient of strain between the notches. Ifju^^ developed a special strain gauge that overcomes many 
of the problems associated with the ASTM recommended strain gauges. Micro-Measurements offers this 
strain gauge configuration as N2A-08-C032A.^^ It is a two-gauge rosette with the gauges at 45° and 
-45° and records the average strain between the notches. These strain gauges were recommended in the 
Micro-Measurements catalog as being specifically made for losipescu shear testing (fig. 30).^^ The 
strain gauge outputs the shear strain directly, avoiding the need for correction factors. The strain gauge 
installation was the same as used for the tensile test specimens. 
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Figure 30. Shear specimen with strain gauges. 


The gauges are mounted on each side of the test specimen in order to remove the effects 
of twisting from the specimen. The strain gauges were connected into a full-bridge circuit as recom- 
mended.^*’23 configuration of four strain gauges, two on each side of the specimen, results in a 
strain reading that is twice the desired reading. A metallic “pathfinder” specimen was used to verify the 
losipescu test fixture and instrumentation (fig. 31). The pathfinder specimen had a known shear modulus 
value of 10.8 x 10^ Ibf/in.^ The percentage of error is shown in the following calculation: 


Error% = x 1 00% = -1 .2% . 

10.8 

The test data compared very well to the actual value for the shear modulus. 


(20) 



Figure 31. losipescu pathfinder specimen. 
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3.7 Tensile Test Procedure 


The tensile test was run using an MTS 20-kip test machine with a numerical control system 
(fig. 32). The strain gauge and extensometer data were recorded with a digital data acquisition system. 
The load was applied at a rate of 0.05 in. per minute and data were recorded every 0.5 sec. Load was 
applied until an auditable “snap” marked the initial failure point of the cementitious material; strain data 
beyond that point were unreliable. 



Figure 32. Tensile test machine. 


As discussed in section 3.5, extensometers and strain gauges were initially used for strain mea- 
surement. The axial extensometer and the strain gauge data were in fairly good agreement (fig. 23). 
However, the extensometers could not be attached to both sides of the test specimens. In addition, the 
knife edges of the transverse extensometers slipped and lost contact with the specimen easily. Therefore, 
the extensometers were replaced with strain gauges. 

Each tensile test specimen had four strain gauges, one axial and one transverse gauge on each 
side. The strain gauges were located back to back on each specimen; this allowed the strain data to be 
corrected for any bending caused by misalignment of the test articles. The data were corrected for 
bending by taking the average of the two strain gauges (fig. 33). 

The strain gauges in the axial direction provided consistent and reliable data throughout the test 
program. However, the strain gauges in the transverse direction proved to be less trustworthy. The data 
illustrated in figure 34 are typical of the transverse data recorded for most of the tensile specimens. The 
strain data appear to be constant for a number of loads; then the strain suddenly increases. This pattern 
repeats itself, forming a curve that resembles a step function. A linear best-fit line through the data will 
be used to approximate the curve. 
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3.8 Axial Properties From Tensile Test Data 


Numerous tensile test specimens were produced and tested while developing the testing proce- 
dures described in this TM. While this trial-and-error procedure was time consuming and labor intensive, 
it eventually produced an acceptable test specimen and procedure. A final set of tensile test specimens 
was produced using the lessons learned from previous experience. 

The test specimens were each tested using the same procedure. Data from four strain gauges and 
the machine load were recorded. The data acquisition system stored the load and strain data in a spread- 
sheet for analysis and plotting. The stress values are calculated from the area of the test specimen and 
applied load: 


Load Ibf 

(Width)(Thickness) in.^ 


( 21 ) 


As previously discussed, the strain values are the average from the back-to-back gauges: 


e 


axial 


pside 1 . pside 2 
axial axial 


2 


(22) 


The data were then plotted to produce a stress-strain curve (fig. 35) for each specimen. 



Figure 35. Typical stress-strain curve for each specimen. 


31 



Young’s modulus is defined as the slope of the linear portion of the stress-strain curve.^ The 
value can be determined by placing a best-fit line through the data points. The slope of the line is then 
taken as Young’s modulus for the specimen (fig. 36). This procedure was repeated for six test specimens. 
A summary of tensile test data indicates consistent results (fig. 37). Appendix B contains the complete 
set of data plots for the tensile test specimens. 


y=0.388x+ 15.812 



Figure 36. Calculation of Young’s modulus for each specimen. 



Figure 37. Summary of tensile test specimens. 
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Young’s modulus values from the six tensile test specimens are summarized in table 4. The 
“pathfinder” specimen CG-SG4 was made thinner than the other specimens. Several of the thinner 
specimens were damaged during the sanding process, so a thicker mold and specimen were later pro- 
duced. The remaining five specimens were made from the thicker mold and are represented in figure 37. 
The average Young’s modulus value for the five specimens is 4.04 x 10^ Ibf/in.^, with a standard 
deviation of 0.20 x 10^ Ibf/in.^ 


Table 4. Young’s modulus from tensile test. 


Specimen 

Width 

(in.) 

Thickness 

(in.) 

Young’s Modulus 
Test Results 
(Ibf/in.^) 

CG-SG4 

1.004 

0.120 

5.00x105 

Dec_28-4 

1.016 

0.190 

3.88x105 

Dec_28-5 

1.006 

0.196 

3.91 X 105 

Dec_28-6 

1.000 

0.216 

4.07x105 

Dec_28-7 

1.003 

0.210 

3.98x105 

Dec_28-8 

1.000 

0.180 

4.40x105 


3.9 Poisson’s Ratio From Transverse Test Data 

The transverse strain values can be used to determine the Poisson’s ratio for the composite 
material. By definition, Poisson’s ratio is given by the equation Vj 2 = - assuming a uniaxial 
stress.^ The transverse strain data from each of the five test specimens is presented in detail in 
appendix B. As mentioned in section 3.7, the transverse strain data presented here are the average 
of two back-to-back strain gauges (fig. 38). In addition, a best-fit line is drawn through the data points. 
The data were fairly consistent from one test article to another (fig. 39); the same “step function” pattern 
was present in all the test articles (sec. 3.7). 

y = -2.9036x- 2.5466 
900 
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■f 500 

l« 

3 300 
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0 

Figure 38. Typical transverse strain curve. 
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Figure 39. Summary of transverse strain data. 


and 


The Poisson’s ratio calculations were made from the following straight-line approximation: 


(7 — 0.388c -i- 15.812 

(23) 

<y - -2.9036ejj.jjjjg - 2.5466 . 

(24) 


Solving the equations for the strain, 


and 


"axial 


"trans 


((7-15.812) 

0.388 


{a + 2.5466) 


-2.9036 

Substituting into the equation for Poisson’s ratio. 


(25) 

(26) 


then 


p 

y _ trans 

o 

^ axial 


(27) 
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( 28 ) 


(0- + 2,5466) 

_ L -2,9036 
^ r (g-i5,8i2) ~ ■ 

0,388 

Poisson’s ratio can now be calculated for each of the five test specimens (table 5), 


Table 5, Poisson’s ratio. 


Specimen 

Experimental 

Dec_28-4 

0.136 

Dec_28-5 

0.135 

Dec_28-6 

0.142 

Dec_28-7 

0.131 

Dec_28-8 

0.144 


The average Poisson’s ratio is 0,137 with a standard deviation of 0,00532, The small standard 
deviation provides confidence that the data are reasonable, 

3.10 Shear Test Procedure 

ASTM standard D 5379-93^^*^^ documents the standard shear test for determining the shear 
properties of composite materials. This test uses a specimen referred to as the V-notched beam (fig, 40), 
This test is often referred to as the losipescu shear test, in honor of its developer. The test specimen is 
placed in a massive test fixture referred to as an losipescu shear test fixture (fig, 41), The test fixture is 
described in detail in the ASTM standards. 




Figure 40, V-notched shear specimen. 
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Figure 41. losipescu shear test fixture. 


The specimen is inserted into the fixture with the notches aligned with the center of the test 
fixture. The two halves of the test fixture are compressed by the MTS testing machine while monitoring 
load. The relative displacement between the two halves of the fixture loads the notched specimen. Strain 
gauges are placed in the center of the specimen to measure the shear strain response during loading. ^*’23 

3.11 Shear Properties From losipescu Shear Testing 

The cost and availability of the strain gauges used on the test specimens limited the number 
of losipescu shear tests. As a result, only six instrumented test specimens were available. The first test 
specimen was a metallic pathfinder specimen that was used to verify that the test procedure, instrumen- 
tation, and fixture could be used to determine the shear modulus for a well-characterized material 
(sec. 3.5). The results from the pathfinder specimen were in excellent agreement with the known values 
for the material. However, of the remaining five composite specimens, only two specimens were suc- 
cessfully tested. The data from the first composite specimen were not usable due to a problem with the 
data acquisition system. The increment for the load data was not properly set during the test. The strain 
gauges on the third and fifth test articles failed to provide usable data during the test. The gauges appar- 
ently debonded from the test specimens. The debonds were most likely caused by poor surface prepara- 
tion during the strain gauge installation. It should be emphasized that the losipescu shear fixture and test 
specimens performed well during testing. The results from the two successful tests confirm this observa- 
tion (fig. 42). The strain gauge on specimen 4 failed at »3,000 ji\ however, the test article continued to 
take load. The strain gauge on specimen 2 continued to provide strain data until the article failed. 
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Shear 4 Shear 2 



The shear modulus is defined as the slope of the linear portion of the shear stress-shear strain 
curve.^’^^ The value can be determined by placing a best-fit line through the datapoints. The slope of the 
line is then taken as the shear modulus for the specimen (fig. 43). This procedure was repeated for two 
test specimens. A summary of shear test data indicates consistent results (table 6). Appendix B contains 
the complete set of data plots for the shear test specimens. 


y= 0 0719x+ 17 258 



Figure 43. losipescu shear test article 2. 
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Table 6. Shear modulus summary. 


Specimen 

Width 

(in.) 

Thickness 

(in.) 

Shear Moduius 
(ibf/in.2) 

Shear 2 

0.458 

0.180 

7.2x10'* 

Shear 4 

0.457 

0.186 

7.8x104 

Average 

- 

- 

7.5x104 


A best-fit line is drawn through the datapoints to provide a means of accurately calculating the 
shear modulus (fig. 43). The shear modulus is given by the slope of the straight line:^^ 


(305t 


Ibf 


Article 2 


m/ 


■ 89 ^) 


(4,000 /t- 1,000 jU) 




m. 


(29) 


(209 


Ibf 


''Article 4 “ 


m. 


- 91 ^) 

in. 


(2,500/t-l,000/t) 


7.8x10^ 


Ibf 


m. 


(30) 


As shown in table 6, the average shear modulus is 7.5 x 10^ Ibf/in.^, with a standard deviation 
of 0.42 X lO'^ Ibf/in.^ The percentage of possible error is determined by the following: 


and 


7.5 X 104 0.42 x 10^ = 7.92 x 10^ ^ 

in.^ in.^ in.^ 


7.92x10^^-7.5x104^ 

Error% = lJi^ = 5.6% 


7.5x104 


Ibf 


m." 


(31) 


(32) 


3.12 Summary of Composite Material Properties 

The material properties for this graphite-reinforced cementitious composite material were 
successfully determined from standard test methods. It was found that the constituent materials can be 
mixed and placed to form laminated composite plates. The plates were cured at room temperature and 
without a vacuum bag. The composite material can be machined into various shapes using standard 
machining tools such as a router. Strain gauges can be applied to the surface of the composite using 
M-Bond 200 and following the standard methods of installation. The results of the material testing 
were consistent and repeatable for both tensile and shear testing. 


The materials property data can now be used as a basis for comparison to theoretically derived 
values. In addition, these values will be used to develop the analytical tools needed to predict the 
response of multilayered beams made from graphite-reinforced cementitious materials. 
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4. RULE OF MIXTURES FOR A COMPOSITE MATERIAL 


The rule of mixtures is investigated as a method to predict the material properties of graphite- 
reinforced cementitious composites. The equations for the effective material properties are derived using 
the assumptions associated with the rule of mixtures. These material properties are then compared to the 
results of the material properties testing. Several derivations are then made to illustrate that these mate- 
rial properties can be used with the elastic curve equations to determine the deflection of simple beams. 
Deflection results from a finite element model and a static test are compared to the analysis using the 
elastic curve equations. The results illustrate the deflections can be predicted by using the material 
property data from the rule of mixtures method. However, several deficiencies and limitations are 
encountered. The conclusion is made that the rule of mixtures can provide a reasonable value for 
Young’s modulus in the axial direction. However, the theory cannot provide accurate values for 
Poisson’s ratio, the shear modulus, or Young’s modulus transverse to the direction of the fibers and load. 

4.1 Material Properties in the Axial Direction 

The rule of mixtures is a method for determining the effective properties of a composite material 
based on the contributions of the individual constituents.^*^ The properties of the composite are related to 
the quantity of the fibers in the matrix. The fiber volume fraction (Ufiber) the matrix volume fraction 
(^matrix) to measure the amount of fiber or matrix in a given volume of material. E in the direc- 

tion of the fiber, E-^, can be determined by assuming the strains in the fiber direction are the same in the 
fibers as in the matrix; i.e., = £matrix- Since the strain in the fiber is the same as the strain in the 

matrix, the sections normal to the fiber direction remain plane after stressing (fig. 44). 



Figure 44. Loading in axial direction. 
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The strain in direction 1 , the fiber direction, is given by the deflection divided by the original 
length, = AL/L. Applying Hooke’s law and assuming both the fiber and the matrix are linear, then 
and “ ^matrix^l- relationship can be compared to two springs acting in 

parallel (fig. 45). The total load (P) is carried by the load in the fiber (/’fiber) the load in the matrix 

('^matrix)' 


f 


Kf 


“l/l/l/H 




“Vi/i/iH 


Figure 45. Spring in parallel. 


The applied force (P), P = -Ffiber + ^matrix’ cross-sectional area (A). The fiber 

stress (ctfibej.) acts on the cross-sectional area of the fiber (Afjber)- The matrix stress (^matrix) 
on the matrix area (A^^atrix)- The force equation can then be written in terms of the stress and area, 
<7fAi = (7fiberAfiber + <^matrix"^matrix- ^iow by substituting the Hooke’s law equations for the stress terms, 
the force equation may be written in terms of strain and Young’s modulus as £’ieiAj =£’fiber^lAfiber 
+ ^matrix^l^matrix- Simplifying the equation results in an expression for Young’s modulus in terms 
of the ratio of the areas: 


FA FA 
£ _ fiber fiber _j_ matrix matrix 


^fiber ^matrix he defined as Fffbgj- Afjbgj-Mj and Aj^atrix/'^n 


(33) 


By definition, the volume fractions must satisfy the equation Ff^ber + Fjjjatrix = T Thus, 
the expression for the effective Young’s modulus in direction 1 can be written in terms of the volume 
fractions: 


F =V F -\-V F 

^ fiber ^ fiber '^matrix ^matrix 


(34) 


4.2 Material Properties in the Transverse Direction 

The effective property in the matrix direction, perpendicular to the fiber direction, can be 
determined by assuming equivalent stress in the fiber and the matrix (fig. 46). 
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Figure 46. Loading in transverse direction. 


In the transverse direction, the cross-sectional areas are equal and the load is constant across the 
area; therefore, o is equal for the fiber and the matrix. Then the strains in the fiber and matrix can be 
written in terms of the stress, = (72/^iber ^matrix ~ '^2/^matrix- Th® matrix transverse 

strains are not assumed to be equal. Now consider figure 46. The total transverse deflection (Aw) is the 
sum of the fiber deflection (AWf^j^gj.) and matrix deflection (AWj^gjj.jx): Aw = Aw^ji^gj. + Defining 

the strains in terms of the fiber deflection, = AWfj^,gj.^Lfjj,gj. and the deflection of the matrix, Cj^atrix 
= The transverse dimension over which the fiber and matrix strains act can be written in 

terms of Lfjber- ~ ^fiber^ Tjj^atrix “ ’^matrix^- The deflection equation can now be defined in 
terms of strain; e 2 W = ffibg,. (l^fiber^) +^matrix(’'^matrix^)- hlow the w factor vanishes from both sides 
of the equation. Hooke’s law equations can be substituted for the strains, O 2 IE 2 = k^jber '^ 2 /^fiber 
+ '^matrix «^2/^matrix- The StreSS term now vanishes from both sides, 1/^2 = ^fiber/^fiber + ^matrix/^matrix- 
Solving the above equation for the £’2 term gives an equation for the Young’s modulus in the transverse 
direction as 


E 

4.3 


F F 

fiber '^matrix 

^ 1 / F +V F 

matrix fiber ''' fiber ^ matrix 

Determination of Poisson’s Ratio 


(35) 


Poisson’s ratio (v) is defined as the ratio of the transverse strain over the strain in the fiber direc- 
tion when stress in the fiber direction is present and all other stress is zero, v ^2 = hi such a case, 

the following definitions can be made: = Cmatrix/^f ^fiber“ ^iber/^1- Recalling that the deflec- 

tion equation was previously defined as £ 2 W = Sfiber^^fiber^) ^matrix^T^matrix^)’ ^he latter can be divided 
by and simplified to the equation £ 2 /£i = %iber/^f(Rfiber) ^matrix/^l(T^matrix)- Substituting the 
Poisson’s ratio (v^ 2 ) defined above in the equations results in the following expression: 

'^12 “ '^fiber^fiber ^matrixh^matrix • (36) 
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4.4 Determination of the Shear Modulus 


To determine the shear modulus (Gj 2 ), one assumes that shear stress (t) on the fiber 
and the matrix are equal. As illustrated in figure 47, the shear for the composite plate is defined 
7 i2 “ ^1^12- 


*■ 


n 


Matrix 

Fiber 

i 

VJ 






Direction 2 

Matrix T 

1 ^ Direction 1 





^ L ► 

Figure 47. Shear loading. 


This assumption means that the shear deformation for the matrix and the fiber can be written as 
/matrix ~ '^/^matrix /fiber = "^/^fiber- shown in figure 48, the shear deflection (A) is a function of the 
shear strain and the width (w): A= yw. In addition, A is the sum of the fiber deflection plus the matrix 

deflection; A = Afiber + A^atrix- 



Figure 48. Shear deformations. 
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As previously shown, the deflection for the matrix and fiber can be stated in terms of the volume 

fractions, and l^matrix' ^matrix “ %iatrix(^matrix^) '^fiber “ Itiber^^fiber^)' Substituting into 
the deflection equations and allowing the tv factor to vanish results in the equation for shear strain: 
y~ Tmatrix^matrix + Tfiber^ fiber- Finally, one can substitute the Hooke’s law equations for shear strain 
to obtain t/G ^2 = Fjj^atrix'^'^'^matrix This equation can be simplified by allowing the shear 

stress term to vanish and solving for G^^ 2 . 


Gi2 = 


c c 

^matrix^fiber 


Y . Gc x. -\-Vr x. G 

matrix fiber fiber matrix 


(37) 


4.5 Calculation of the Effective Material Properties 

Equations (34)-(37) can now be used to calculate the effective material properties of a graphite- 
reinforced cementitious composite. The material properties for the matrix were obtained from the 
manufacturer of the graphite mesh. The properties for the matrix were obtained from previous research^ 
and confirmed with tensile tests (app. B). The material properties for the fiber and the matrix are defined 
as follows: 


^fiber 

= 33.5x106 lbf/in.2 

p 

matrix 

= 115.1 xl03 lbf/in.2 

'^fiber 

= 0.3 

'Vatrix 

= 0.27 

Pfiber 

= 0.064 lbf/in.2 

Pmatrix 

= 0.35 lbf/in.3 

pfiber 

= •^fiber/2(l + Vf^ber) 


^matrix ^matrix(^(^ '^matrix^- 

The fiber area in the composite must be determined from the properties of the graphite. 

Data provided with the graphite define the area in one graphite yam The number of yams 

(Ay am) is determined for the test specimen and the total area of the graphite (A^^j^gj.) is determined: 

Ayarn = 0.000175 in.2, Ayam = 8, and Afi^er = iVyam'4yarn- 

The total area of the specimen (Aj) is determined from the thickness and width of the specimen. 
For the example presented here, the test article is labeled CG-SG4, and t = 0.12 in., iv = 1.004 in., and 
Ai= 0.012 in.2 -pj^g matrix volumes, respectively, can be calculated as follows: Tjjjjgj. = 

and Tfiber = 0 . 012 ; ^matrix = 1 - Ffiber and Vmatrix = 0 - 988 . 

Thus, the values for the effective material properties can be calculated as follows: 
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Young’s modulus (direction 1): 


F =V F -\-V F 

^ fiber ^ fiber ^ matrix ^ matrix 


=5.045 xl05lbf/in.2 


Young’s modulus (direction 2): 


E2 = 


F F 
^ fiber ^matrix 


V F -\-V F 
matrix ^ fiber ^ fiber ^matrix 


£2 = 1 . 164x105 lbf/in.2 


Poisson’s ratio: 


^12 ^fiber^fiber '^matrix^matrix 


V 12 - 0.27 


( 38 ) 


(39) 


(40) 


• Shear modulus: 


^12 


C C 
^fiber^matrix 


4^matrix^fiber ■*" ^^fiber^ 


matrix 


G ^2 = 4.62x104 lbf/in.2 


(41) 


4.6 Comparison of Rule of Mixtures Properties to Test Data 

Prior to each tensile test, the width and thickness of each specimen was recorded. Specimen 
CG-SG_4 was the pathfinder specimen and was made thinner than the other specimens. The hand- 
polishing of the specimen surface resulted in a variation in specimen thickness for the remaining five 
specimens (table 7). These dimensions were used with the Mathcad® solution sheet for the rule of 
mixtures to predict the Young’s modulus in the axial direction. The results of the analysis are presented 
in table 7. The results from the rule of mixtures can be compared to the results from tensile testing. 
Figure 49 is a stress-strain curve for specimen CG-SG_4. 
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Table 7. Young’s modulus predicted from the rule of mixtures. 


Specimen 

Width 

(in.) 

Thickness 

(in.) 

Young’s Moduius 
Ruie of Mixtures 
(ibf/in.^) 

CG-SG4 

1.004 

0.120 

5.05 xIO® 

Dec_28-4 

1.016 

0.190 

3.58x105 

Dec_28-5 

1.006 

0.196 

3.52x105 

Dec_28-6 

1.000 

0.216 

3.31 X 105 

Dec_28-7 

1.003 

0.210 

3.37x105 

Dec_28-8 

1.000 

0.180 

3.75x105 



Figure 49. Composite stress-strain curve (SG_4). 


The testing of this type of specimen is discussed in detail in section 3. The data are presented 
here to provide a comparison between the material properties calculated from the rule of mixtures and 
the test data. Young’s modulus from the tensile test is calculated by taking the slope of the stress-strain 
curve: 
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(1,500-500) 


^test “ 


Ibf 

in.2 


,5 Ibf 


(3,000- 1000) X 10 




^test - 5 X 10 2 

m. 


^error 


^test ^1 


"test 


xlOO% 


Aerror = 0.898% . 


( 42 ) 


The analysis is repeated for each of the tensile test specimens. Comparing the Young’s modulus 
data to predicted values provides insight into the accuracy of the rule of mixtures method (table 8). 


Table 8. Comparison of predicted Young’s modulus with test results. 


Specimen 

Young’s Modulus 
Rule of Mixtures 
(lbf/in.2) 

Young’s Modulus 
Test Results 
(lbf/in.2) 

Measurement of Error 

(%) 

CG-SG4 

5.05x105 

5.00x105 

0.89 

Dec_28-4 

3.58x105 

3.88x105 

7.70 

Dec_28-5 

3.52x105 

3.91 X 105 

10.00 

Dec_28-6 

3.31 X 105 

4.07x105 

18.60 

Dec_28-7 

3.37x105 

3.98x105 

15.30 

Dec_28-8 

3.75x105 

4.40x105 

14.80 


The errors range from 0.89 to 18.6 percent; averaging the six test specimens yields an average 
error of 1 1.23 percent with a standard deviation of 6.4 percent. In each case, the rule of mixtures under- 
estimates the test value, suggesting that the predictions may be modified^’^*^ to more closely match the 
test values. Empirical modifications to the rule of mixtures would require a considerable test program to 
generate a large enough database from which to derive an accurate factor. It is sufficient for this TM to 
recognize that the rule of mixtures provides a reasonable estimate of the composite material properties. 
The comparison establishes the fact that the rule of mixtures is applicable to a cementitious -based 
composite and suggests that employing empirical modifications can increase the accuracy. 

Young’s modulus in direction 2 was not accurately determined from the rule of mixtures. This 
composite is made from a graphite mesh, which results in fibers running in both directions 1 and 2. 

Since the composite has the same material composition in both directions, it is intuitively obvious that 
the Young’s modulus in directions 1 and 2 would be equal. The rule of mixtures calculates the E 2 value 
well below test results. This is not unexpected, since this theory considers only fibers in direction 1.^40 
However, the equations developed by KrencheP can be used to transform the fiber direction and load to 
accurately predict £’ 2 - KrencheTs transformation equations depend on the orientation of the fibers to the 
load. These equations are based on the same principles as presented in section 4.2. However, the proper- 
ties in direction 2 can be calculated by mathematically rotating the load 90° (app. A). Regardless of the 
apparent inadequate results in direction 2, the rule of mixtures can still provide the necessary data for the 
analysis of simple beams. 
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4.7 Comparison of Poisson’s Ratio From Tensile Test Data 

Poisson’s ratio from the transverse data is presented in section 3.8. Poisson’s ratio for the five test 
specimens presented in table 4 is calculated using the rule of mixtures equation (40), v ^2 = ^fiber^fiber 
+ %atrix^matrix- "*"^0 results of the calculations are presented in table 9. 

Table 9. Comparison of predicted Poisson’s with test results. 


Specimen 

Experimental 

Rule of Mixtures 

Error 

(%) 

Dec_28-4 

■■ 



Dec_28-5 



48.5 

Dec_28-6 

0.142 

0.262 

45.8 

Dec_28-7 

0.131 

0.262 


Dec_28-8 

0.144 

0.262 

■■ 


The error in the Poisson’s ratio calculation is substantial. However, it is not entirely unexpected 
(sec. 3.5).^’^*^ The derivation for the rule of mixtures assumes the reinforcement fibers are present in the 
axial direction only. The theory cannot account for fibers running in the transverse direction, emphasiz- 
ing the need to develop a solution based on the laminated plate theory for composites. 

4.8 Comparison of Shear Modulus Values 

The shear modulus (sec. 4.4) of the matrix and the fiber can be obtained knowing the Young’s 
modulus and Poisson’s ratio as follows:^ 


G 


matrix 


E 

matrix 

2(1 + '^matrix) 


4.5 XlO^^ 
in.^ 


Tiber 


^fiber 


2(l + Vfiber) 


= 12.8x106 


Ibf 


m." 
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The shear modulus of the composite is determined from the shear moduli of the constituents 
and the volume fractions of the composite. Applying equation (41), 


Gi2 = 


c c 

^matrix ^fiber 


V C -\-V c 
^ matrix'^fiber fiber '^matrix 


• = 4.6x104 


Ibf 


m." 


(44) 


It is apparent from the G ^2 value that the graphite fiber does not significantly increase the shear modulus 
of this composite. However, the results of the shear testing do not support this observation (fig. 40). 
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Comparing the predicted value and the value from the test article (table 10), 


(4.6x104-7.2x104) 

Error% = ^ ^ x 100% = 36. 1 % . 
7.2x104 


( 45 ) 


Table 10. Comparison of predicted shear modulus with test results. 



Shear Modulus 
Rule of Mixtures 

Shear Modulus 
Test Results 

Measurement of Error 

Specimen 

(lbf/in.2) 

(lbf/in.2) 

(%) 

Shear 2 

4.6x10^ 

7.2x104 

-36.1 

Shear 4 

4.6x104 

7.8x104 

-41.0 


The average percent of error, 38.6 percent, indicates that the rule of mixtures significantly 
underpredicts the shear modulus for this composite. The test results indicated the shear modulus is 
greater for the composite than for the cementitious material. 

4.9 Summary of Material Properties 

In the axial direction, Young’s modulus can be determined from the rule of mixtures. The rule 
of mixtures prediction is on average 1 1 percent less than the value obtained from tensile test data. As 
expected, the rule of mixtures did not predict the Young’s modulus in the transverse direction. For the 
composite considered here, F’axiai ~ ^transverse- ^ addition, the rule of mixtures did not predict the 
Poisson’s ratio for the composite material. The material property testing indicated a Poisson’s ratio value 
=50 percent less than the value predicted by the rule of mixtures. The results of the shear test determined 
a shear modulus that is on average 39 percent higher than the value predicted by the mle of mixtures. 
Although only two shear test specimens were available, the good correlation in the two test specimens 
gives credibility to the test data. 

The limitations of the rule of mixtures are obvious from the review of these results. This observa- 
tion leads to the conclusion that additional analytical methods must be mvestigated to determine all the 
material properties for the composite. However, before pursing these additional analytical methods, it 
should be recognized that many analytical problems could be solved using the information provided by 
the rule of mixtures. 

4.10 Calculation of the Deflection for a Beam in Pure Bending 

The rule of mixtures method for calculating the effective material properties for a composite 
beam has been compared to material test data. The results indicated that the rule of mixtures could be 
used to predict the Young’s modulus in the direction of the fiber and load, or F^effective- 
logical step is to combine this material property with an analytical method that can be used to predict the 
behavior of a composite structure. As noted earlier, the transform section method was not accurate for 
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this composite because the Young’s modulus ratio was extremely high. However, using the effective 
materials property data, the Young’s modulus ratio (rj = ^’effective/^matrix) <20 and in the accept- 

able range. The Biszick beam^ can now be analyzed by using the effective material property data and the 
transform section theory. The results will be compared to static test results to give an indication of the 
accuracy of the methods employed. 

The cross section of a composite beam, fabricated of the same graphite and matrix materials used 
throughout this research, is illustrated in figure 50. An assumption is made that the beam is composed of 
two materials; the “graphite section” is composed of graphite and matrix material. The hg dimension of 
this section is the same thickness as the tensile test articles used to verify the effective material proper- 
ties for the composite material. The transform section analysis method will assume the graphite- 
reinforced section has the effective material properties calculated by the mle of mixture method. The 
geometry of the test specimen is used to determine the dimensions of the graphite strands and the thick- 
ness of the graphite section: 

= 1 .20 in. 
h = 0.115 in. 

o 


Graphite Section 


Matrix Section 


Figure 50. Cross section of beam. 



The remaining dimensions of the composite section are calculated as follows: 

= 0.361 in. - hg 
w^= 1.20 in. 

= 2.46 X 10-1 in. 

This test specimen has nine yams across the width. The fiber volume can be determined 
from the geometry of the test specimen: 

Ayj^rn = 0.000175 in.2 
N =9 

yarn 

Afu =N A 

nber ’ yarn yarn 

A^ =(0.115 in.) (1.20 in.) 

^fiber “ ^fiber/^1 

Yfibe,= 1.141x10-2. 
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The matrix volume is then calculated from the rule of mixtures: 


^matrix ^ ^fiber 
^matrix = 9-886 X lQ-1. 


The material properties of the graphite and matrix are the same as previously defined: 
= 33.5 X 10^ Ibf/in .2 


^matrix ' 


115x103 lbf/in.2 


The effective material properties for the composite section of the beam can be calculated. 
Young’s modulus in direction 1: 


F =V F 4-V F 
1 fiber ^ fiber matrix matrix 


= 4.9602x103 7 


Ibf 


(46) 


mF 


The moment of inertia for the beam can now be calculated using the transform section method. 
The Young’s modulus ratio for the transform section is calculated from the effective material properties 
and the matrix material properties: 


V = 


77 = 4.313 . 


(47) 


' matrix 


The beam geometry is transformed by multiplying the width of the graphite section by rj: 


w = rjw 


g 


w = 5.176 in. 

6 


(48) 


The centroid and the area of the graphite section are given by 


y<, = K + 


fhA 


V ^ / 




(49) 


while the centroid and the area of the matrix material is 


3'r=- 


4 = 


(50) 
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The centroid of the transform section is given by 




=2.4366 xlO-lin. 


(51) 


The cross section of the composite is symmetric in the X direction; thus, the X location 
of the centroid is the center of the beam; 


w 

y = 

^bar 2 


^bar = 0-6 in. 


(52) 


The moment of inertia for the graphite-reinforced section of the beam can be calculated by using 
the parallel axis theorem, 1 = 1+ A (fi. The base of the top section is transformed by rj and the height of 

o o 

the graphite-reinforced section is given by h = h . The moment of inertia of the graphite-reinforced 

o 

section with respect to its centroid is given by the equation 


_h/;3 
\2 ' 

The area of the graphite-reinforced section is the base multiplied by the height: 

A ^=bl, . 

The distance from the center top section to the centroid of the section is 


(53) 


(54) 




d = Y, 


bar 


h 




V ^ y 


(55) 


Applying the parallel axis theorem. 


total 


's_, 01.1=2-7875x10-3^4 


(56) 


The moment of inertia for the matrix (bottom) section of the beam can be calculated using 
the same method. The base and the height of the matrix section are 


b = w^ 


h = hc . 


(57) 
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The moment of inertia of the matrix section with respect to its centroid and the areas are 


bh^ 


12 


A^ = bh . 


( 58 ) 


The distance from the center of the matrix section to the centroid of the section is 


h \ 
^ 7 


(59) 


The moment of inertia for the matrix with respect to the centroid of the composite beam 
can be calculated as follows: 


total 


‘c ,„B,= 5^7864x10-3111^“ 


(60) 


The total moment of inertia for the transform section is the sum of the moment of inertia for each 


section: 


^total _ total _ total 


= 8.5739 xl0-3in.4 


(61) 


The beam was subjected to a three-section bending test. A rigorous deflection analysis of the 
beam test was performed to establish a base of comparison to finite element modeling and test data. 
The comparisons provide insight into the accuracy of the methods used for the analysis. 


The equation for the elastic curve of this beam will be derived as part of this analysis. Although 
a solution to this beam problem is available in the literature, it was important to develop the solution in 
this TM because of variations in the boundary conditions. Laboratory testing strives to duplicate the 
exact boundary conditions; however, it is sometimes impractical. The equations aided in determining 
which boundary conditions had significant effects on the results being recorded. These equations were 
solved for a variety of boundary conditions using a Mathcad^^ solution sheet (app. C). It should be noted 
that with the aid of these equations and WD-40® lubricant, the test fixture eventually duplicated the 
proper boundary conditions. 

In general, the elastic equation for a curve is given by El((f-ldx^)y = M(x). The moment in each 
section of the beam (fig. 51) can be derived by passing a section through that location. In section 1 
(x < a), for example, the moment at x is computed from the free-body diagram shown in figure 52 as Px. 
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Figure 51. Free-body diagram of the beam. 



Figure 52. Free-body diagram of section 1. 



(62) 


(63) 


(64) 
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Similarly, for section 2, figure 53,x> a but x<b: 



p U 

Figure 53. Free-body diagram of section 2. 


Finally, for section 3, figure 54, j > (a + b), 

r/2 

EI^y. = 2Pa + Pb-Px 
dx^ 

d Px^ 

El—y. = 2Pax + Pbx + Cc 

dx 2 ^ 



(65) 

( 66 ) 

(67) 


(68) 

(69) 

(70) 
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Summarizing the equations and replacing the derivative with 6, 


Px 

EIQiyi = -- + Ci , 


(71) 


EI62y2 = + C3 , 


Px^ 

ElO^y^ = ^Pax + Pbx — — + , 


(72) 

(73) 


Px^ 

Ely I = + C^x + C 2 , 


(74) 


EIy2 = —ax^+C2X + C^ , 


(75) 


and 


Phx'^ P 

Ely 2 = Pax^ + — — + C 5 X + Cg 


(76) 


The equations must satisfy the boundary conditions. For example, at x = 0, = 0, equation (74) 

can be written 


0 = -x^+C^Xo+C2 . 


(77) 


At x = a, equation (74) must equal equation (75); therefore. 


P . P . 

x^ + E^x^ + C2 — dXi + C'lX^ + C4 . 
6 2 


(78) 


In addition, at v = n, = y 2 , and = 02 - Equation (71) must equal equation (72): 


Px^ 

— — + Cj = Paxi + C3 


(79) 
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For the boundary condition x = (a + b), equation (75) equals equation (76): 


— ax2 + C^X2 + 6 * 4 = Pax^ + ■ 


Pbx'^ Px} 


+ C^X2 + Cf^ . 


( 80 ) 


Continuing in the same manner, at x = (a + b), J 2 = ^ 3 , and 02 = ^ 3 » then equation (72) equals 
equation (73): 


Paxn + Co = 2Pax^ + Pbx^ 


Pxi 


■ + C. 


(81) 


Finally, for the boundary condition x = L, equation (76) must equal zero: 


« Pbx% Px\ 

0 = Pax^ H J- + C5X3 + C^ 


2 


(82) 


Equations (71)-(76), with the given boundary conditions, represent a system of equations which 
can be solved simultaneously for the constants Cj through Cg. A Mathcad® solution was developed to 
perform the calculations (app. C). Appendix C contains typical values for these constants. The elastic 
curve equations may be restated as a function of j: 

In section \,x<a: 


and 




Px^ 


+ C^x + C2 


e^(x) = 


^ Px^ 


El 


+ Ci 


(83) 


(84) 


The value for the moment of inertia is set to the value calculated from the rule of mixtures and 
the transform section method. The Young’s modulus in the equation is set equal to that of the matrix; i.e.. 


1 = 1 


total 


/ = 8.574 xl0-3in.4 


E = E 


matrix 


C = 1.15x105 


Ibf 


(85) 


m.^ 
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In section 2, x> a < a + b: 


and 


Finally, in section 3, x> a + b: 


y2(x) = ^(Pax^+C^x + C4) 




( 86 ) 


(87) 


/ 




Pax^ + 


V 


Pbx^ Px^ 


+ C^x + Cg 


and 


( 88 ) 


03(^) = 


El 


Px^ 

IP ax + Pbx 1- C;: 

2 ^ 


(89) 


The composite beam was then subjected to the three-section bending test (fig. 55). Deflection 
data were recorded at the center of the 9-in. beam. 



Figure 55. Three-section bending test fixture. 
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The elastic beam analysis was repeated for each of the load steps used in the three-section beam 
test. The load and deflection results are listed in the following matrices: 


0 


0 


0 

3.43 


-0.0445 


-0.0420 

5.64 


-0.0729 


-0.0690 

7.41 


-0.0940 


-0.0910 

9.62 

Ibf, Jtest = 

-0.1197 

in., }’cal = 

-0.1180 

11.84 


-0.1420 


-0.1450 

14.05 


-0.1670 


-0.1720 

16.26 


-0.1920 


-0.1990 

17.81 


-0.2090 


-0.2180 


(90) 


Figure 56 illustrates that the test results and the analytical data are in very good agreement. Error 
bars on the analytical data represent the 6.4-percent standard deviation associated in determining 
Young’s modulus (table 8). This indicates that the transform section method of analysis can be used to 
predict the deflections. However, the rule of mixtures must be applied to determine the effective Young’s 
modulus. 


— Test Data — B— Analytical Data 
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4.11 Development of a Finite Element Model 


An additional point of comparison was obtained by developing a finite element model of the 
three-section beam test (fig. 57). The finite element model was developed for two reasons: (1) To con- 
firm that the elastic beam analysis was correct and (2) to gain insight into modeling composite structures 
using finite element codes. The finite element model was executed for each of the static test loads. 



z 



X 


Default Deformation: 
Max2.46-0f@Nd B40 


Figure 57. Finite element model. 


The deflection results from the finite element model “specimen 1” are as follows: 


Jfem 


0 

-0.0420 

-0.0690 

-0.0910 

-0.1180 

-0.1450 

-0.1720 

-0.1990 

-0.2180 


m. 


(91) 
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The deflection results for both the elastic beam equation and the finite element model were plotted 
with the deflections from the three-section bending test. 


It is significant that the transform section method of analysis returns the same deflections as the 
finite element solution. This leads to the conclusion that the transform section method of analysis can be 
used for the analysis of a composite in which the Young’s modulus ratio is high, provided the effective 
material property data are included in the calculations. Figure 57 indicates that both of the analytical 
methods satisfactorily predict the results of the bending test (fig. 58). The error between the elastic beam 
equations and the test data is relatively small: 


^ >'test step->'2(4-5in.) 

terror = X 100% 

•f test _ step 


A 


error 


-4.52% . 


(92) 


— Test Data —B— Analytical Data A Finite Element 



Figure 58. Transform section, finite element, and test results. 


4.12 Calculate the Deflection for a Five-Layer Plate in Pure Bending 

Section 4.2 demonstrates the analysis of a composite beam made from the graphite -reinforced 
cementitious material. The University of Alabama in Huntsville (UAH) concrete canoe team produced 
a large plate that is similar to the layups used in canoe construction. A test article 5 in. wide and 15 in. 
long was cut from the composite plate. The challenge is to determine the deflections of the plate using 
the information gained from the material properties of the constituents, the mle and mixtures, and the 
transform section method of analysis. These are the tools and information that are available to the 
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students at the time they are designing these competitive canoes. The beam is a buildup of five layers 
as shown in figure 59. 



Figure 59. Cross section of multilayered composite beam. 


The properties of the beam can be calculated by using the transform section method and the 
effective material properties as determined from the rule mixtures. The effective Young’s modulus for the 
three graphite and matrix sections is 


IK'f 

£ =5.04x105— r 
^ in.^ 


(93) 


The Young’s modulus for the matrix material, F'j, is determined from tensile testing of the 
cementitious material (app. B); 


^115x103-^ (94) 

in.^ 

It was previously established that the Young’s modulus ratio must be <20 for the transform 
section method to be accurate. Using the values established here, the Young’s modulus ratio is calculated 
and is within acceptable limits: 



77 = 4.383 . 


(95) 
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The dimension of the graphite sections is transformed by multiplying the width of the sections by 
rj. A Mathcad solution sheet (app. C) was used to perform the calculations for determining the moment 
of inertia for the beam. These calculations are also described in section 4.10. The moment of inertia for 
this cross section was calculated to be 


= 38.2x10-3 in.4 (96) 

In general, the procedure for determining the properties of multilayer beams is the same as that 
for a single-layer beam. It is interesting to note that the graphite section in the middle of the beam does 
not contribute significantly to the bending strength. The beam is placed in a three-section bending test 
with the specimen supported over a 15-in. span and loaded at two points, each located at a distance of 
2.25 in. from the center of the span. Members of the UAH canoe team performed the testing during the 
1999-2000 school year. This test was typical of those used by the team to evaluate the composite 
designs considered for the canoe. 

The elastic beam equation for this problem was provided previously. The deflection of the center 
of the beam was calculated for the same loads used in the beam tests. Equations (86) and (87) are used 
for the deflection calculations: 


y2ix) = -^{Pax^ + C2X + C4) ( 9 ' 7 ) 

and 

e2ix) = ^(Pax + C^) . ( 98 ) 

The calculations were made using the Mathcad solution sheet (app. C). The derivation 
for the beam equation is described in detail in section 4.10. Young’s modulus used in the equation 
is from the matrix material: £’^ = 115x103 Ibf/in.^ 

The results of this analysis are compared to the results from the bending test. As shown in figure 
60, they compare well at deflections <0.25 in., then diverge as the deflections reach higher values. The 
test results also appear to be slightly nonlinear. The nonlinear behavior of the displacements suggests 
that the material properties are nonlinear. 

The stress-strain curve for the test specimen can be established from the flexure formula 
and the moment in the beam. The strain is calculated from the stress and based on the Young’s modulus 
of the matrix as follows: 


(J(x) = 


M{x)c 

I 



e{x) = 


o{x) 

F 

matrix 


(99) 
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Deflection (in.) 


Figure 60. Deflection of multilayer beam versus test data (linear). 


These calculations are presented in detail in the Mathcad solution sheet (app. C). The stress- 
strain curve from the test specimen can now be compared to the stress-strain curve for the matrix 
material. Figure 61 shows the stress-strain curve for the beam and the matrix material. 



Microstrain (ju in./in.) 


Figure 61. Stress-strain curve for matrix material. 
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The experimental curve in figure 61 is derived from a combination of test and analysis. The first 
section of the curve is from tensile testing using extensometers in the axial direction; the tensile test 
specimens failed at ~2,000 microstrain. The stress-strain curve beyond that point is derived from the 
tensile test data using a least-squares method to determine an equation for the curve.^ 


n n . \2 ^ 

-t ^ ^matrix ; ^2 S (‘^matrix J “ S ^matrix ^ 

i=l i=l i=l 



These equations are solved for the minimum values of hg, bp and bp 

- Z>o + ^i^^niatrix ^2 (^matrix) * (101) 

Equation (101) is the resulting stress-strain relationship used to extend the data for the tensile 
test. The resulting curve is an approximation that does not represent actual material properties and is 
only used to facilitate the analysis. This analytical method was developed by Balaguru^ for the nonlinear 
analysis of ferrocement in bending. 

The question arises as to whether the tension or compression properties of the matrix are deter- 
mining the nonlinear behavior of the plate. Therefore, the compression properties of the matrix must also 
be considered.^ Figure 62 compares the analytical curve for tensile properties with the compression 
properties from the compression test. The data for the tensile (analytical) and compression (test) proper- 
ties are in good agreement up to 6,000 microstrain. Furthermore, both curves indicate nonlinear behavior 
beyond the 2,000 microstrain level. 

Since Young’s modulus is in the nonlinear region, it must be recalculated as the load increases 
beyond 2,000 microstrain. The elastic beam equation is then used with the nonlinear Young’s modulus 
to calculate the deflections at the higher load levels. These calculations are illustrated in the Mathcad 
solution sheet (app. C). The results are compared to the test data from the three-section bending test. The 
error bars on the analytical data represent the 6.4-percent standard deviation associated in determining 
Young’s modulus. As illustrated in figure 63, the data compares fairly well to the test results, even in the 
nonlinear region of the curve. 
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Figure 62. Tension versus compression properties for matrix material. 
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Figure 63. Comparison of test results with nonlinear transform section results. 






The preceding analysis indicates that the combination of effective material properties and the 
transform section method can be used to determine the deflection of multilayered composite beams. This 
is similar to the work presented by Balaguru^ for the nonlinear behavior of ferrocement. These data also 
indicate that as the matrix material becomes nonlinear, the deflection curve follows the same trend, 
indicating that the bond between the fiber material and the matrix material may be slipping. However, 
it should be noted that the beam did not fail even when the tension strain levels were far above the 
maximum tension strain found when testing the matrix material. 

It must be pointed out that several deficiencies exist in the preceding analysis. The first is recog- 
nizing that the Young’s modulus ratio is constantly changing as the matrix material becomes nonlinear. 
This implies that the moment of inertia for the beam is also changing. These changes are not considered 
in the analysis. In addition, this method of analysis does not provide a means to determine the strains or 
stresses in each layer of the composite beam. 

4.13 Summary of the Rule of Mixtures Analysis 

The rule of mixtures provides a method of calculating the effective material properties of the 
composite from the volume fraction and material properties of the constituents. Comparing the effective 
material properties with the data from the tensile and shear tests provided a measure of accuracy for this 
method. Young’s modulus in the axial direction (£'j) was 1 1 percent lower than the value found in the 
tensile test. However, Young’s modulus in the transverse direction (£’2) and Poisson’s ratio (v^2) did not 
agree with the test data. The poor comparison for £2 ^d Vj2 was not unexpected, since the composite 
under consideration has fibers running in the transverse direction. These fibers are not considered in the 
rule of mixtures theory. The shear modulus Gj2 was 39 percent lower than the value found in two test 
specimens. 

Although the information derived from the rule of mixtures is limited, it was used successfully to 
predict the deflection of single and multilayer composite beams. The analysis used the transform section 
method and the elastic beam equation to demonstrate that the effective Young’s modulus (£^) could be 
used to predict deflections for a beam in pure bending. The analysis was verified by comparison to test 
data and finite element models, and the analysis was extended to predict the nonlinear behavior of 
composite material. 

The fact remains, however, that these derivations are limited to the case in which the fibers are 
running in the axial direction. This limitation can be overcome by applying the laminated plate theory 
of composites to the cementitious composite in question. 
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5. LAMINATED PLATE THEORY OE COMPOSITES 


The rule of mixtures is limited by the assumption that the fibers are aligned in the direction of the 
load Laminated plate theory for composites offers a more complete analytical tool, provided this 
class of material is compatible with the assumptions and derivations involved with the theory. Beginning 
with laminated plate theory, the laminated composite plate equations will be derived. The derivations 
will be programmed into a Mathcad^"^ solution sheet (app. D) to provide an automated method for 
solving the equations. The equations will be arranged to form a system of simultaneous equations from 
which an exact solution can be found. The solution includes all elements of the stiffness matrix and can 
be used for symmetric or nonsymmetric laminated plates. The program is used to analyze a composite 
plate for which the results are well known. Comparisons are made with other programs to demonstrate 
an increase in accuracy. A finite element model is generated to provide a verification of the methods and 
establish the procedures for modeling composite beams. 

5.1 Plate Theory 

Figure 64 shows a load-free plate where the midplane point A coincides with the x-y plane. The 
deflection in the z direction is zero for all points on the plate. The deflection components at a point in the 
X, y, and z directions are given by u, v, and w, respectively. 



Figure 64. Load-free plate. 


The fundamental assumptions of the small deflection theory of bending for isotropic, homoge- 
neous, elastic thin plates are based on their geometry. These assumptions are known as the Kirchhoff 
hypotheses and can be stated as follows:^® 

• The midsurface deflection is small compared with the plate thickness. The slope 
of the deflected surface is very small and the square of slope is negligible. 

• The midplane remains unstrained subsequent to bending. 
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• Plane sections initially normal to the midsurface remain plane and normal to surface 
after the bending. 

• The stress normal to the midplane (o^) is small compared to the other stress components 
and may be neglected. 

A plate can be deformed by normal, shear, and bending loads. The displacements in the jc, y, z 
directions are defined as u, v, w, respectfully. Figure 65 shows a plate experiencing normal stresses.^® 
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Figure 65. Plate with normal displacements. 


The strains for this displacement state in the x and y directions are defined by 


and 



(102) 



(103) 


respectively. 

Figure 66, on the other hand, shows that when shear stresses are applied to a plate, they produce 
strains that are a function of both u and v displacements.^® 

The equation for the in-plane shear strain is 


xy 


d d 

— u-\ V . 

dy dx 


( 104 ) 
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Figure 66. Shear displacement of plate. 


The total in-plane displacement at any point in a plate is the sum of the normal displacements plus the 
displacements induced by bending. Figure 67 shows the geometry associated with a plate in bending.^ 





Figure 67. Bending displacement. 


The displacement along the z direction due to bending is z x sin(0). Since 6 is small, sin(0) = 6, 
and the displacement is z6. In addition, from the geometry, it is clear that 6 is the slope of the plate, 
didx X w. 
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The displacements of the midplane of the plate in the x and y directions are defined as Uq and Vq, 
respectively. Then the displacement in the x direction can be defined in terms of the midplane displace- 
ment, the slope of the plate, and the location from the midplane as 

d 

U=U()-Z—Wr,. (105) 

dx 

Similarly, the displacement in the y direction is 

d 

= (106) 

The Kirchoff hypotheses state that plane sections remain plane. This implies that there is no 
strain in the z direction. According to Kirchoff, the normal and shear strains in the z direction may be 
written as 


£ y = W £7=0 

A 1 A 

dz 


d d 


rxz = o 


d d 

dy dz 


ryz = ^ 


(107) 

(108) 
(109) 


When the displacements in equations (105) and (106) are substituted into equations (102)-(104), 
the strain equation accounts for the midplane displacements and the bending displacements and is 
written as 






( 110 ) 


J2 


= Vn - Z 

y dy ^ dy^ ° 


( 111 ) 


rxy=— “0+— ^0-2z 


dy ^ dx 


d d 
dx dy 


-Wr 


( 112 ) 
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These equations can be further simplified by defining the midplane strains as 


_ d 

“0 

d 



0 


(113) 

(114) 


T^Oxy 


— Un + • 


dy 


dx 


(115) 


Substituting equations (1 13)-(1 15) into equations (1 10)-(1 12) results in the strain equations written 
in terms of the midplane strains 


d r/2 

^ dx ^ dx^ ° 

r/2 


(116) 


(117) 


r/2 


y dy ^ dy^ ° 
d^ 


(118) 

(119) 


d d ^ d d 


( 120 ) 


Yxy roxy ^^dxdy"^^ 


( 121 ) 


The midplane surface curvatures are now defined as follows: 


.10 




0 


( 122 ) 




K^, = 7rWn 

y ^1,2 0 


dy" 


(123) 
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(124) 


i. ^ d d 
K^, = 2 Wrt 

^ dxdy ^ 


Finally, the strain equations may be stated in matrix form: 






f 

K 

1 


= 


— z 


^xy 


^Oxy 


1 


(125) 


5.2 Generalized Hooke’s Law for Nonisotropic Materials 

For isotropic materials (properties are the same in any direction), the relationship between stress 
and strain is independent of the direction of the force. Young’s modulus is the only elastic constant 
required to describe the stress-strain relationship for a uniaxial force. 

For a nonisotropic material, at least two elastic constants are needed to describe the stress-strain 
behavior of the material. Consider, for example, the case of a nonisotropic plate made of a fiber- 
reinforced material (fig. 68). The fibers are in direction 1, resulting in the direction 1 properties being 
much higher than the direction 2 properties. 



Figure 68. Nonisotropic plate with fiber reinforcement. 


A special case of anisotropy is the existence of two perpendicular planes of symmetry in material 
properties. Such materials are referred to as orthotropic. This configuration is characteristic of a lami- 
nated composite that may have different material properties in the perpendicular planes. In general, the 
properties of the material are direction specific. 


72 



The analysis of an orthotropic system must take into account the direction of the material proper- 
ties. For such a system, the stress can be expressed for each direction as 

CTj = (126) 

and 

C72 = £’ 2£2 • ( 127 ) 

Consider the case shown in figure 69; an applied load acts parallel to the fibers. 



Figure 69. Applied load parallel to fibers. 


This type of loading produces stress in more than one direction within the plate and is referred 
to as plane stress. The plate is elongated in direction 1 and contracted in direction 2. This illustrates the 
Poisson’s effect, which may be defined as 


Vl2 = 


z£2 

ei 


(128) 
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The total strain in the orthotropic material must consider the Poisson’s ratio effect for each 
direction. The strain in the plate is now a function of several material properties: 


and 



^ 12^1 


( 129 ) 



^ 21^2 • 


(130) 


Shear forces must also be considered in the plate analysis. Shear stress is related to the shear 
strain by a material constant called the shear modulus (G). As with the normal stress, the shear stress ( t) 
must be defined in terms of the direction of the material as follows: 


^12 “ ^ 12^12 * 


(131) 


In general, the equations for strains can be expressed in terms of stress and the orthotropic 
material properties (Young’s modulus, Poisson ratio, and the shear modulus) as 


2 


( 7 , (T 


(132) 


P -^2 ^ 

^2 -T; '^12'^ 

^2 ^1 


(133) 


ri 2 


T 2 

^12 


(134) 


These equations can be expressed in matrix form as follows: 


1 

1 



_L _ni 0 " 
£, E2 

_L 0 

£i E2 

(^1 

Lri2j 


0 0 

Jl2_ 



Gi2 



( 135 ) 
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The 3x3 matrix in equation (135) is made up entirely of material properties for the orthotropic material. 
This matrix is called the compliance matrix (5): 


1 


El 



0 


Y2L 

El 

J_ 

0 


0 

0 

1 

G\i 


Equation (135) may be rewritten to form an equation for the stress: 


(136) 


^2 

7 

II 

f 

1 

_^12_ 


bi2j 


(137) 


Furthermore, Vi2/£i = V 21 IE 2 , then the stiffness matrix {Q) can now be defined as the inverse of the 
compliance matrix (5) as follows: 


Q = [S]-^ = 


-Vl2^1 


(-•£’2 +^12^2 i £' i ) 
El 


{-E2 + ViiViiEi 


-yiiEi 


El 


{-E2+Vi2V2\Ei 

-1 


(-E2+V12V21E1) 


(138) 


G^12 


Finally, the stress equation can be expressed as a function of the stiffness matrix (Q) to further simplify 
the equation: 


^1 



^2 

= Q 

^1 

_^ 12 _ 


Jl2_ 


(139) 
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5.3 Plane Stress for Generally Orthotropic Plates 


The previous discussion was limited to the case in which the load is either parallel or perpendicu- 
lar to the principal material directions of the lamina; i.e., reinforcement fibers. In general, the loading is 
not in the principal material direction. The equations for the generally orthotropic plates transform the 
stress and strain into coordinates that coincide with the principal material directions. Figure 70 is a free- 
body diagram of a composite laminate; it indicates the directions of the stress in relation to the direction 
of the reinforcement fiber. Direction 1 is the direction of the fiber and direction 2 is perpendicular to the 
fiber. 



Figure 70. Direction of stress. 


The transformation of plane stress is a fundamental part of mechanics of materials. In the case 
of composite materials, the idea is to transform the stress from the direction of the load to the principal 
material direction. As shown in figure 71, the element is cut to produce an area, AA. A free-body 
diagram is used to sum the forces in direction 1 . 



Figure 7 1 . Sum force in direction 1 . 


76 


The sum of the forces in direction 1, for example, results in the following equation for cr^: 


0 = c^AA - f7^(AAcos0)cos0- fJ^(AAsin0)sin0 

+ -T^(AAcos0)sin0-T^(AAsin0)cos0 . (140) 

From the free-body diagram of figure 72, the sum of the forces in direction 2 results 
in an equation for '• 


0 = (72 AA - (7^ (AA sin 0) sin 0 - (7^ (AA cos 0) cos 0 
+ T^(AAcos0)sin0 + T^(AAsin0)cos0 . 


(141) 



Figure 72. Sum force in direction 2. 


Using figure 72, sum the forces in direction 1 to obtain an equation for t^ 2 - 

0 = Ti 2 ^ + cr^AAsin0cos0- ^CT^AAcos0sin0j 

+ |-T^AAcos0cos0j + T^AAsin0sin0 . (142) 

Simplifying equations (140)-(142) and solving for the stresses in the direction of the principal material 
direction yields the stress transformation equations: 

= (7^cos0^+(7^cos0^+2T^sin0cos0 (143) 

O 2 = sin02 + Oy cosO^ - 2t^ sin0cos0 (144) 

t^2 = • (145) 
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Equations (140)-(142) can now be written in matrix form: 


1 

1 


^2 

= 

1 

to 

1 



COS0^ 

sin0^ 

-sin0cos0 


sin0^ 
cos 6^ 
sin0cos0 


2sin0cos0 
-2sin0cos0 
cos0^ -sin0^ 



The 3x3 matrix is known as the transformation matrix (T): 


( 146 ) 


T = 


cos 6^ 
sin 6^ 

-sin0cos0 


sin0^ 2sin0cos0 

cos0^ -2sin0cos0 

sin 0COS 6 cos 9^ - sin 9^ 


Thus, the stress transformation equation ( 146) may be written in terms of (T) as follows: 


(147) 


^1 



^2 

= T 


1 

i 




(148) 


Similarly, to transform from the 1-2 coordinate system to the x-y coordinate system, use the inverse 
of the transformation matrix: 


1 

Q 

1 


1 

1 


II 

to 

1 

i 


\jn\ 


(149) 


The strains can be transformed in the same manner as the stresses. Equations (150) and (151) represent 
the transformation of strain: 





^2 

= T 

s 

ri2 


y xy 

_ 2 _ 


_ 2 _ 


( 150 ) 
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( 151 ) 


■ 







= r-i 

^2 

'yxy 


ri 2 

_ 2 _ 


_ 2 _ 


However, the strain transformation equation (151) creates a numerical problem. The shear strain 
term ( 712 /^) compatible with the shear strain term from equation (139), Yi 2 - To overcome this 
problem, a factor of 2 is introduced into the matrix equations as follows: 


- 







1 0 0 



s 

= 

0 1 0 



1 


0 0 2 


'^xy 
_ 2 _ 


(152) 


The resulting matrix is called Reuter’s matrix (R): 


R = 


1 0 
0 1 
0 0 


0 

0 . 
2 


(153) 


The shear strain values can be multiplied as needed to combine the transformation equations (148)- 
(151) with the compliance equation (139). These derivations are provided in the Mathcad program listed 
in appendix D. The resulting equation is as follows: 


<y^ 



X 


X 


= T-^QRTR-^ 


7xy_ 


'y xy 


(154) 


Equation (154) illustrates the transformed stiffness relationship between the stresses and strains 
for an orthotropic plate. This equation provides the stress-strain relationship in the x-y direction. 
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Equation (154) is difficult to work with because of the number of matrices involved. However, 
the equation can be simplified to three matrices by defining the transformed reduced stiffness matrix 

(Q): 


Q =T-^QRTR-^ 


(155) 


Finally, equation (154) may be written in terms of Q : 






= Q 




'^xy 


(156) 


The Q matrix is a 3 x 3 matrix that represents the transformation of the stiffness matrix (Q). The terms 
in the matrix contain four independent material properties and the angle between the reinforcement 
fibers and the load. 


5.4 Mechanics of Laminated Composite Plates 

The strain equations for a plate in bending must now be related to a plate made from several 
layers or lamina. Solving equation (156) for the strain values results in the following: 





X 

= Q-^ 

X 



^xy 


_^xy_ 


(157) 


The midplane strains and the curvatures can now be equated to the stress values by combining 
equations (125) and (157), 


1 

Q 
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1 

O 
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1 
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^Oy 

— z 


f 
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^Oxy 


>1 
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(158) 
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Simplifying equation (158) and solving for the stress matrix results in an equation for the stress in terms 
of the midplane strains and curvatures: 






K, 


= Q 

^Oy 

-zQ 

Ky 

'^xy 


YOxy 


Kycy 


(159) 


Equation (159) represents the stress in each ply of the laminate for an orthotropic plate. 

5.5 Determination of Force and Moment Resultants 

The resultant forces and moments acting on a laminate are a function of the stresses on each 
layer through the thickness of the plate. The force and moment resultants are important for the analysis 
of composite plates because the stress in the plate varies from top to bottom. The force and the moment 
act at the midplane of the plate in the same direction as the stresses. 

Figure 73 illustrates the normal and shear forces that act on a plate. The resultant force can be 
found by the summation of the stress multiplied by the area. The sum of the forces can also be found 
by integration of the stress multiplied over the thickness (t) of the plate, then multiplied by the area. 
Figure 74 represents a cross section of a plate. The width (w) times an incremental thickness (dz) 
represents the area over which the force is acting. 



Figure 73. Resultant forces. 



Figure 74. Integration over width of plate. 
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Equation (160) represents the forces acting on the plate: 


t 

^ o^w dz = dz . 

« T 

The width (w) of the plate occurs on both sides of the equation and cancels. This allows 
the resultant force to be defined in units of force per unit length. The resultant force per unit length 
at the midplane of the plate is defined as 


= (IM) 

2 

Similarly, for the remaining directions of stress. 


and 



In matrix form, these equations can be written as 


(162) 


(163) 



t 


X 


X 




y 

J-t 

3 ' 


2 

1 

1 


(164) 


Figure 75 illustrates the moments acting on a plate. The resultant moments for the midplane 
of each plate can also be calculated from the stresses. The moments are defined as the sum of the stress 
multiplied by the area, multiplied by the moment arm with respect to the midplane. The moments are 
also defined in terms of unit lengths. The moment equations are as follows: 
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(165) 


dz 

2 
t 

OyZ dz (166) 

~2 

t_ 

2 



Figure 75. Resultant moments on plate. 


In the same manner as the forces, the moments can be stated in a matrix equation: 


t 




r — 
2 


My 

= 



M^, 


-r 

2 


xy 

J 

L 


5.6 Laminate Constitutive Equations 


(168) 


Equations (164) and (168) define the forces and moments acting on the midplane of the plate 
in terms of the plate stresses. Equation (159) defines the stresses acting on any layer in the plate in terms 
of the midplane strains and plate curvatures. The midplane strains are functions of the midplane dis- 
placements and the midplane curvatures are functions of the displacement (w). Equations (159), (164), 
and (168) are the plate constitutive equations. 
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Equations (164) and (168) must be modified to define the midplane forces and moments for each 
layer in the plate. This can be done by recognizing that the stress resultants can be expressed as the sum 
of the stress resultant for each layer. Letting n represent the number of layers in the plate, the equation 
for the forces becomes 




'h 


n 



=i 




k=i 






xy 

J 

h{k-l) 

L TVj 


(169) 


Similarly, the equations for moments are written as 


r -| 
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xy 
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h{k-\) 

xy 


(170) 


The integration parameters have now changed from the thickness of the plate to the location 
of the top and bottom of each layer. Figure 76 illustrates the integration terms. 



Figure 76. Layer location. 
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The equation for the transformed stress, equation (159), can be substituted into the equations for 
the midplane forces and moments. The resulting equations (171) and (172) state the forces and moments 
in terms of the strain and curvatures: 





r- 

p -| 


p 


n 



^Ox 


K, 

^y 

II 


Q 

^Oy 

+ zQ 


Nxy_ 

j 

h{k-l) 

_ 

^Oxy 


Kxy_ 


dz 


(171) 
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(172) 


Integrating these equations is relatively simple because the matrices £q and K are not functions 
of z and are not included in the integration. The plate stiffness matrix ( 2 ) is constant over the thickness 
of any layer; equation (172) can be written as 
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Performing the integration results in the following equation: 
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(173) 


(174) 


85 



The moment equation can also be simplified as 
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( 175 ) 
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4-4-1 


(176) 


Note that the midsurface strains (£q) and the curvatures (K) are not part of the summations. 

The stiffness ( Q ) and the hj^ terms can be summed outside equations (174) and (176) to form 
new matrices, and the equations become 


and 
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(177) 


(178) 


In equation (177), matrix A is called the extensional stiffness matrix. This matrix relates the 
normal stress and strains similar to the modulus of elasticity. However, if the matrix terms A^g and A 23 
are nonzero and the plate has a shear load applied to it, normal stresses will be produced: 


A = 


n 


Qk,Mk~^k-l) ■ 


( 179 ) 


The thickness of the ^h layer can be defined as where tj^ = hj^~ ^ Then equation (179) can be 
written as a function of the thickness of the layer: 
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(180) 


k=\ 


The B matrix in equations (177) and (178) is called the coupling stiffness matrix. The B matrix 
relates bending strains with normal stresses: 




k=\ 


hk-hl-i 


(181) 


This equation can also be stated as a function of the layer thickness (z^) and the distance from the layer 
midplane to the center of the ^h layer. Defining z^ with the following equation, 


_ h+h-\ 

The equation for the coupling matrix B may be stated as 


(182) 


n 

^ = X Qk,\ h ^k • 


k=\ 


(183) 


The coupling stiffness matrix has substantial influence on a laminated composite plate. A non- 
zero B matrix implies that an axial load will cause the plate to bend and/or twist. Also, any moment 
applied to the plate will cause extension of the plate midplane. This coupling between bending and 
extension is key to understanding composite plates, and neglecting the effects of the B matrix could be 
catastrophic to a composite structure. 

The D matrix in equation ( 178) is called the bending stiffness matrix. This matrix relates 
the plate curvatures with the bending moments: 


k=\ 

The equations for D can also be written in terms of and z^ : 


o=Set,i 
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(184) 
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5.7 Effective Engineering Properties for Laminated Composite Plates 

The effective engineering properties are determined from the A, B, and D matrices. The concept 
is to determine the longitudinal Young’s modulus (EJ, transverse Young’s modulus (Ey), shear modulus 
G^y, and longitudinal Poisson ratio (v^) of the laminate (plate) from the E 2 , G 12 , and v ^2 material 
properties of one of the layers. A plate can be designed and optimized by varying the angles of the 
reinforcement fibers in each layer. 

Recalling Hooke’s law for an isotropic plate, E^ = oje^. Similarly, the relationship developed 
for laminates relates the midplane force to the axial strain: 


E 


X 



(186) 


where h is the laminate thickness. 

Equations (177) and (178) may be combined to form a linear system of equations. The force 
and moment equation may be written as follows: 
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(187) 


This relationship is a system of six equations and six unknowns. The midplane strains and curvatures 
are the unknowns. The A, B, and D matrices are constants. 

The force and moment matrix can be defined to determine the effective material properties 
for the composite. In the symmetric case, the layers above the midplane are identical to the layers below 
the midplane, resulting in zero values for the B matrix. This greatly reduces the size and complexity 
of the solution. For this analysis, however, the case of nonsymmetric plates is considered. 

5.8 Calculation of Effective Material Properties 

The solution to the system of equations (187) is best demonstrated through an example prob- 
lem. A Mathcad solution sheet (app. D) contains the calculations used in this example. Consider, 
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for example, a plate that is made up of four fiber-reinforced layers as shown in figure 77. The angles 
of the reinforcements are 0° and 45° and each layer is 0.005 in. thick. 



Figure 77. Angle of fiber reinforcement. 


The following material constants are assumed for an arbitrary composite material: 

• Young’s modulus in the principal material direction 1: = 20.01 x 10^ Ibf/in.^ 

• Young’s modulus in the principal material direction 2: £’2 = 1.301 x 10^ Ibf/in.^ 

• Shear modulus in the principal material direction 1-2: Gj 2 = 1.001 x 10^ Ibf/in.^ 

• Poisson’s ratio in direction 1-2: v ^2 = 0.30 

• Poisson’s ratio in direction 2-1: V 21 = v^ 2 ^ 2 l^\- 

The compliance matrix is determined as follows from equation (137): 

0 

0 
1 

The stiffness matrix, inverse of the compliance matrix, is calculated from equation (138): 


4.998x10-8 -1.499x10-8 0x10^ 

-1.499x10-8 7.686x10-^ 0x10^ 


OxlOO 9.99x10-7 
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2.013x10^ 3.926x105 
2 = 5-1= 3.926x105 1.309x106 
0 0 


( 189 ) 


1.001x106 


The geometry of the layup is then used to determine the transformation matrix. These data are presented 
in detail in the Mathcad solution sheet (app. D). The transformation matrix is derived from the stress 
transformation equation. The matrix is solved for each ply angle: 


cos 02 
sin 02 


sin02 2sin0„ cos On 
cos 02 -2 sin On cos 0^ 


- sin 0^ cos On sin 0^ cos 0^ cos 02 - sin 02 


1 0 0 


0.5 0.5 1 


Plyl:Ti= 0 1 0 Ply 2: 72= 0.5 0.5 -1 

0 0 1 -0.5 0.5 0 


(190) 


0.5 0.5 1 


1 0 0 


Ply3:Ti= 0.5 0.5 -1 Ply4:Ti= 0 10. 

-0.5 0.5 0 0 0 1 


(191) 


As previously described, Reuter’s matrix (R) is used to condition the Q matrix: 

"l 0 0“ 

/?= 0 1 0 . 

0 0 2 


(192) 


The Q matrix represents the stiffness and transform matrix; the Q matrix is calculated for each 
ply using equation (193): 

Q=T-^QRTnR-^ ■ ( 193 ) 


2.01x102 3.93 xl05 

0 

3.93x105 1.31x106 

0 - 

1 0 

0 

lxl06 

6.56x106 4.55x105 

4.7x106 

4.55x105 6.56x106 

4.7x106 

4.7x106 

4.7x106 

5.16x106 


(194) 


( 195 ) 
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Ibf 


(196) 


2s 


6.56x106 4.55x106 4.7xl06 

4.55 xl05 6.56x106 4.7xl06 

4.7x106 4.7x106 5.16x106 



2.01x10'^ 
3.93 xlO^ 
OxlOO 


3.93 xl05 
1.31x106 
0x10° 


OxlOO 

OxlOO 

lxl06 



As expected, the plies with the same angles are equivalent. 


The A matrix, or the extensional stiffness matrix, is calculated using equation (180): 


(197) 
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k=i 


2.67 xlO^ 
4.95 xlO^ 
4.7 xlO^ 


4.95 xl04 
7.87 xl04 
4.7x10^ 


4.7x10^ 

4.7x10^ 

6.16x10^ 


Ibf 


in. 


(198) 


The B matrix, or the coupling stiffness matrix, is calculated using equation (183). In the case 
of a symmetric plate, the B matrix is zero: 
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Qk^k^k 
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0-0 0 
-0 -0 0 
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M 

in. 


(199) 


The D matrix, or the bending stiffness matrix, is calculated using equation (185). These values 
indicate how the bending moments will affect the plate curvatures: 
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1 2, 
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IL+t 7^ 
19 ^h^k 
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6.09x10-1 1.31xl00 
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3.92x10-1 
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( 200 ) 


Since the A, B, and D matrices are defined, equation (187) can be solved for the strains 
and curvatures. However, the loads must first be defined. In the Mathcad solution sheet, the load 
matrix is created to represent the forces and moments at the midplane of the laminate. The load matrix 
is from equation (187). In the first case, a uniaxial load is applied in the x direction: 
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(201) 


A global matrix that represents the stiffness matrix in equation (187) is constmcted from the A, 
B, and D matrices. Equation (187) may now be written in terms of the loads and the global stiffness 
matrix: 
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(202) 


This system of equations represents a plate with a uniaxial load applied in the x direction. 

There are six equations and six unknown strain and curvature values. A Mathcad subroutine, Isolve, 
is used to find a solution to the system of equations. The Isolve function returns values for strains 
and curvatures. Hooke’s law can then be applied to determine Young’s modulus and Poisson’s ratio: 

• Young’s modulus in jc: 

= -^ = 1.13x10^^ (203) 

• Poisson’s ratio: 

V =-^=3.17x10-1 . (204) 

The effective modulus in the y direction is found by assuming a unit force in the y direction and 
setting all other forces and moments to zero. The system of equations is now solved with the boundary 
conditions for a uniaxial load in the y direction. Then based on Hooke’s Law, Young’s modulus and 
Poisson’s ratio can be calculated: 
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Young’s modulus in x: 


= 2.0972 X 106 ^ (205) 

- h£oy 

• Poisson’s ratio: 

=^^ = 5.87x10-2 . (206) 

The effective shear modulus is found by assuming a unit force in the x-y direction and setting all 
other forces and moments to zero. The shear modulus (G„,) is a function of the shear load and the shear 
strain: 

• Shear modulus: 




^>r0xy 


3.22 X 106 

in.^ 


(207) 


5.9 Comparison of Results for the Symmetric Solution 

The engineering properties for a symmetric laminate can be compared to solutions from a num- 
ber of sources. This comparison illustrates that the method used for finding a solution to the constitutive 
equation (187) significantly affects the accuracy of the calculations. For example, Halpin^ reduces the 
constitutive equations by assuming the Aj 3 and A 2 3 terms in the A matrix are zero. This results in the 
following equations for the material properties: 
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X __ Hatpin 





^2,2^^ 


(208) 


and 


"y_ Hatpin 


A.t'^2,2 (^t,2) 

Atj/l 


(209) 


G 


xy _ Hatpin 



( 210 ) 


Nettles, on the other hand, includes all the terms in the A matrix and is able to develop the following 
equations for the materials properties: 
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^2,3 


h 


~^2,3 ^ A, 3 A, 2^3, 3 (^,3) ^2,3 
'^ 3,3 21^^(213 3) -(^13) 


( 212 ) 


and 


G 


xy _ Nettles 


^ 3,3 

(^2,3) , 2^13^1,2^2,2^2,3 (^1,2) 

(^2,3) (^,3) (^2,2) 

h 

^'^^ 2,2 /, 

A,l(^ 2 , 2 ) 

(^1,2) ^ 2,2 



( 213 ) 


Although Nettles’ method is an improvement over the Halpin method, these equations do not apply 
to the nonsymmetric case. When the 3 and A 2 3 terms are set equal to zero, the Nettles^® equations 
reduce to the Halpin^ equations. 

Table 1 1 is a comparison of results obtained from several commercially available computer 
programs for laminated composite plate analysis. 


Table 11. Results from different laminated plate solution techniques. 



Ex 


^xy 

Exact 

1,13x10^ 

2.1 x1Q6 

3.22x106 

Nettles 

1,13x10^ 

2.1 x1Q6 

3.22x106 

Halpin 

1,18x10^ 

3.47 x1Q6 

3.08x106 

Inplane 

1,18x10^ 

3.48 x1Q6 

3.08x106 

Patran 

1,18x10^ 

3.48 x1Q6 

3.08x106 


As expected, the Nettles equations provide the same solution as the exact solution. However, 
these equations do not work for the nonsymmetric case. The Halpin, Inplane, and Patran solutions 
overestimate the stiffness of the composite in the transverse direction and underestimate the torsional 
stiffness for the composite. 

It is difficult to determine the specific effects this type of inaccuracy could produce. It is the 
intent of this study to indicate that these inaccuracies can be avoided by using the solution techniques 
developed in the previous sections. 
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5.10 Determination of Strain 


Section 5.9 determined the effective material properties of the laminated composite plate. It is 
also desirable to determine the stress and strain in each layer of the plate. This is accomplished by using 
constitutive equations (177) and (178) to determine the midplane strains and curvatures. The total strain 
in each ply or layer is then found by using equation (125). The resulting strains are then transformed by 
equation (150) to the principal material direction. As shown previously, the constitutive equation can be 
solved for the midplane strains and curvatures. 

For example, assume the laminated plate in the first example is loaded with 1,000 Ibf in the 
X direction. The plate is 5 in. wide. The x direction is parallel to the 0° plies. The calculations for this 
example are contained in the Mathcad solution sheet (app. D). 

As stated in the problem definition, a force is applied in the x direction. The plate is 5 in. wide 
so would be calculated as follows: 




1,000 Ibf 
5 in. 


The load matrix in equation (187) will contain the following values: 


(214) 




’ 2 x 102 ‘ 



OxlOO 



OxlOO 

^x 


OxlOO 

My 


OxlOO 



0x10® 


(215) 


The midplane strains and curvatures can be determined in the same manner as demonstrated 
in section 5.9. The solution indicates that the midplane strain includes axial, transverse, and shear 
strains, while curvatures are zero: 

Midplane Strains: Midplane Curvatures: 

£(3^ = 8.82x10-4 /«:^=0xl00 

£(3^ =-2.8x10-4 /«:^=0xl00 

yo„=^.6xio-' Ar^=0xi0». 

Equation (125) can be used to transform the midplane strains and curvatures into the strains in each ply. 
To determine the strains in the first 45° ply, define Zpjy as the distance from the midplane to the center 
of the ply: z-^^y = 2.5 x 10-^ in. 
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The strains in the first 45° ply in the x-y coordinate system are calculated from the midplane 
strain and curvatures (eq. (125)). The values are then transformed into the principal material directions 
for the 45° ply with equation (150). The details of the calculation can be found in the Mathcad solution 
sheet (app. D). The resulting matrix is the strain in the principal material direction: 




7.13x10-5 

^2 

= 

5.31x10-4 

/12_ 


-1.16x10-5 


5.11 Determination of Stress 


(216) 


The stress values can be calculated from the strain by applying equation (155) to determine 
the Q matrix for the given ply. The Q value for the 45° ply is as follows: 



6.56 x1Q6 
4.55 x 1Q6 
4.7x10^ 


4.55 x1Q6 

6.56 x 1Q6 
4.7 xlO^ 


4.7x10^ 

4.7x10^ 

5.16x10^ 



(217) 


The strain in the 45° ply for layer 2 is determined from equation (125). Equation (156) is then used 
to determine the stresses in the ply: 




«x 


2.35x105" 




= 

1.99xl0i 



y xy 


4.6x102 



Now use equation (148) to transform the ply stress into the principle material direction: 


(218) 


^1 




1.64x105 " 

^2 

n 


= 

7.23x102 

_'^12_ 


'^xy 


-1.16x105 



(219) 


5.12 Determination of Displacement 

The deflection of the plate is found from the midplane strains and the curvatures of the plates. 
Since these values are known, the deflections can be found by integration: 
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( 220 ) 


Ax= r 

J y 




Jr, I J r 


-K^dx 


\dx . 


The limits of integration are along the edge of the plate. In the example, the plate is symmetric. As a 
result, the coupling stiffness matrix and the curvatures are zero. The Mathcad solution sheet (app. D) 
includes the calculations; the curvature values are included for the nonsymmetric case. In the example 
problem, the deflection in the x direction is calculated as 

Ar = 8.82xl0-^in. (221) 

5.13 Determination of Neutral Axis 

The neutral axis is the plane of zero strain for any direction in an unsymmetrical plate. 

The location of the plane may be determined from equations (125) and (187). For pure bending 
in the x direction, the force and moment matrix is zero except for the moment about the a: axis: 


0 

0 

0 

Mx 

0 

0 


^11 

^12 

^13 

^11 

^12 

D \3 

^21 

^22 

^23 

^21 

^22 

^23 

^31 

^32 

^33 

^31 

^32 

^33 

^11 

^12 

^13 

Du 

A 2 

^13 

^21 

^22 

^23 

^21 

D22 

D23 

.^31 

^32 

^33 

D31 

D32 

D33 


^Ox 

^Oy 

^Oxy 

Kx 

Ky 


(222) 


Solving for the midplane strains and curvatures and then reducing the matrix equation yields 
the following expression: 


X 

0 

1 


1 

1 

^Oy 


^21 

'^Oxy 


^31 

Kx 


^11 



^21 

1 


_^ 31 _ 




(223) 


The above equation can then be separated to form an expression for the midplane strains and curvatures. 
Equation (125) can then be used to find the location of the neutral axis. These derivations are contained 
in the Mathcad solution sheet (app. D). 
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5.14 Comparison of Results to a Finite Element Model 

A finite element model of a simple plate is developed to verify the previous set of equations. 

The plate is 10 in. long and 5 in. wide. As in the previous example, the plate has four layers oriented 0°, 
45°, 45°, and 0°. As shown in figure 78, the plate is clamped at the left end with a load of 1,000 Ibf 
applied at the right. The model uses MSC/Nastran as the finite element code. 



Figure 78. Boundary conditions on finite element model. 


The input for the composite material properties was calculated using equation (187). The effec- 
tive material properties from equation (187) for the model were defined in the MSC/Nastran format by 
using a MAT2 card. The input for the MAT2 card is represented by the following matrix equation: 


^2 

= 

1 

Gi2 

^22 ^23 


^^1 

^^2 

1 

1 

«2 

_'^12_ 


L^31 

^32 ^33 _ 


Jl2_ 


«3 


(224) 


For the purpose of this example, the thermal effects are neglected. The above equation becomes 
the stress-strain relationship of equation (139). As previously discussed, Q is the stiffness matrix where 


Q = 


j_ 

El 

-V21 

El 

0 


-V21 

E2 

J_ 

El 

0 



0 

1 


Gii J 


(225) 
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The material property data for the finite element code are calculated using equation (225). 

The MAT2 input card is written to reflect the following Q: 

Gi 2 G13] [2.013x10'^ 3.926x105 0 1 

Q= G 21 G22 G23 = 3.926 x1Q5 1.309x106 0 ^ (226) 

G31 G32 G33J 0 0 1.001x106 

The properties of the composite plate are defined by the PCOMP command. This input is unique 
to the MSC/Nastran code and requires no calculations. The input is listed for completeness as follows: 

PCOMP, 1, , , 10000.0, STRN, 

+, 1, .005, 0.0, YES, 1, .005, 45.0, YES, 

+, 1, .005, 45.0, YES, 1, .005, 0.0, YES . 

The ply stress and strain in the principal material directions are calculated using the Mathcad 
solution sheet with the geometry and loading from the example problem. The results are provided and 
compared to the results of the finite element model below. 

These results are the ply strains in the principal material direction for the 45° ply: 

• Laminated plate theory results: 

"gj 1 r 7.08x10-5 " 

£2 = 5.283x10-4 (227) 

712J -1.166x10-5 

• Finite element model results: 

"ci] [7.13x10-5" 

£2 = 5.31x10-4 . (228) 

712J -1.16x10-5 

These results are the ply stress in the principal material direction for the 45° ply: 

• Laminated plate theory results: 

"(Ti] [ 1.633x105 

(J2 = 7.194x102 (229) 

T12J -1.1673x105 


99 



Finite element model results: 


^1 


"i.64x103 

^2 

= 

7.23 xl02 

Jl2_ 


-1.16x10^ 


(230) 


The deflection of the plate in the direction of the load (x) compared very well with the finite 
element model: 


Laminated plate theory results: 




Ax = 8.8x10-3 

(231) 

Finite element model results: 




Ax = 8.8x10-3 . 

(232) 


The comparison of the results shows that the finite element model agrees perfectly with the Mathcad 
solution for the plate under uniaxial load. 

5.15 Conclusions From Laminated Plate Theory 

As stated previously, the purpose of the laminated plate derivations and the development of the 
Mathcad and finite element solutions is to create methods for the analysis of complex multilayered, 
graphite-reinforced cementitious composite structures. The development of the Mathcad solution sheet 
proved that there is no numerical constraint to the laminated plate theory as found in the transform 
section theory (sec. 1.4). In addition, it has been determined that the material property data developed 
previously are adequate for the analysis of any conceivable combination of layers and graphite orienta- 
tions. An improvement has also been made to the accuracy of the solution method used to solve the 
laminated plate equations. 
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6. ANALYSIS OF MULTILAYERED SYMMETRIC AND NONSYMMETRIC 
GRAPHITE-REINFORCED CEMENTITIOUS BEAMS 


Section 6 demonstrates that cementitious composite material can be analyzed using the laminated 
composite plate theory. Section 5 verified the laminated plate equations for an arbitrary composite by 
comparison to the results of a finite element model. The equations will now be applied to the analysis 
of three graphite-reinforced cementitious beams subjected to pure bending. 

The rule of mixtures is used to determine the effective and E 2 values. The shear modulus 
and Poisson’s ratio values are taken from the results of the material properties tests. The laminated plate 
theory is then used to determine the effective material properties for each of the beams. These properties 
are used with the elastic beam equation to predict the deflection of each beam. The results from the 
analysis are compared to the results from the test. In addition, the laminated plate equations are used 
to derive a set of material constants for a finite element model. The results from the model are also 
compared to the test results. 


6.1 Introduction 

As shown in figure 79, each beam will consist of two graphite-reinforced layers. The graphite 
mesh in the first beam will not be rotated. In the second beam, both layers of the graphite mesh are 
rotated 45°. In the third beam, the bottom layer of graphite is rotated 45° while the top layer will not 
be rotated, producing a nonsymmetric laminated section. 



Figure 79. Cross section of multilayered composite beam. 
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In each case, the laminated plate equations are used to determine the effective material constants. 
Elastic beam equations, developed in section 4.1 1, are used to determine the deflection at the center of 
each beam. The analysis is repeated for several load increments to produce a load versus deflection plot. 

A finite element model is developed for each beam consistent with the methodology developed 
in section 4.12. The MSC/Nastran finite element program is given the data for each layer in the form of 
MAT2 and PCOMP input commands. The finite element model computes the deflection at the center of 
the beam, and the model is run with incrementing loads to produce the load versus deflection plot. 

Finally, a comparison is made between the analytical data and the test data, and conclusions 
are drawn based on these results. 

6.2 Three-Section (Pure Bending) Testing 

The beam is placed in a three-section bending test with the specimen supported over an 8.75-in. 
span and loaded at two points, each located 1.375 in. from the center of the span. This is the same three- 
section bending test that is described in section 4.10. 

The deflections were measured at the center of each beam using a digital deflection gauge 
(fig. 80). A cup was hung from the loading hook of the fixture. Loading was achieved by incrementally 
placing a measured amount of lead pellets in the cup. The time between load increments was not con- 
trolled but could be estimated at =30 sec. The displacement was recorded at each load increment. Each 
beam was tested by the same method and with the same loading increment. 



Figure 80. Three-section bending test. 
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The results of the beam tests, illustrated in figure 81, show that beam 1 [90,90] is much stiffer 
than beam 2 [45,45]. As expected, beam 3 [45,90] has a stiffness value between beams 1 and 2. The 45° 
layer is on the tension side of the pure bending test, which may explain why the beam behavior is close 
to beam 2. 


Beam 1 [90,90] » Beam 2 [45,45] a- Beam 3 [90,45] 

1,600 
1,400 
1,200 
3 1,000 

V) 

(/) 

O) 

53 800 

1 600 

03 

m 

400 
200 
0 

0 0.02 0.04 0.06 0.08 0.1 0.12 

Deflection (in.) 

Figure 81. Deflection plot of multilayered beams. 



















































6.3 Calculation of the Deflection Using Elastic Beam Equations 

The deflection of the center of the beam was calculated for the same loads used in the three- 
section bending tests. The elastic beam equation for this problem was provided in section 4.10. 

The moment of inertia is calculated by / = M^/12. The Young’s modulus is calculated from the lami- 
nated plate theory. The Mathcad solution sheet (app. C) is used to perform the calculation to determine 
the deflection at each load increment. 

6.4 Determination of Effective Material Properties 

The material properties for each layer of the composite are determined from the rule of mixtures 
and the laminated composite plate theory. For example, using the equations developed in section 4.1 1, 
the material properties of the fiber and matrix are as follows: 


£^,^^^ = 33.5x106 lbf/in.2 

^matrix =115.1x103 lbf/in.2 (233) 
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The area of the graphite material is calculated from the manufacturer’s data as 


Ayarn = 0.000175 in.^ iVyam = 8 

^fiber ~ -^yam '^yarn • (234) 

The area of each layer is the thickness multiplied by the width Ay = 0.120 in. x 1 in.: 


biber=^ = 0012 

= l-b, be, = 0.988 . (235) 

Thus, the values for the effective material properties can be calculated for each layer: 

• Young’s modulus (direction 1): 


p —V F 4-V F 

^ 1 ^ fiber ^ fiber ^ matrix ^ matrix 

Ey = 5.045 X 10^ lbf/in.2 (236) 

This procedure is used to determine the effective Young’s modulus for each layer in the composite beam. 

The shear modulus and Poisson’s ratio values are determined from the material test results 
(sec. 4.14). These effective properties are used with the laminate composite plate equations to determine 
the effective material properties for the composite beam. 

The Mathcad solution sheet (app. D) was set up for a two-ply plate consistent with the procedure 
for determining the effective properties of a multilayer composite plate, described in detail in section 
5.9. The bottom entry in each property array corresponds to the bottom ply of the plate and the geomet- 
ric center of the beam is considered zero: 


5.0x10^ 

Ibf 

5.0x10^ 

Ibf ^ 

7.5 xlO"^ 

Ibf 

0.1 5" 

5.0x10^ 

•2 ^2 - 

m. 

5.0x10^ 

. 2 ^12 “ 
m. 

7.5 xlO"^ 

.2 ^12 = 
m. 

0.15 
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The orientation, thickness, and centroid of each ply are as follows: 


0.125' 


‘ .0625 ■ 

n 

■ 90 " 


in. thickness z = 


in. centroid 0 = 


0.125 


-.0625 

180 

45 


ply angles • 


(238) 


The thickness of each layer is one half the beam thickness. It should be noted the exact thickness of each 
beam varies (0.24 to 0.27 in.) due to imperfections in making the beams. The centroid of each layer is 
measured from the geometric center of the beam. 

The compliance matrix is then calculated for each layer in the composite. The stiffness matrix 
( QJ is then formed from the inverse of the compliance matrix: 


5n = 


1 


El 

-Vl2, 

El 

0 


-^21, , 

Eln 

1 

El 

0 


0 

1 


T2„ 


2n=[S„r‘ . 


(239) 


The subscript n refers to the layer number. 

A stress transformation matrix is calculated from the angles given in the array: 


cos 02 
sin 02 

-sin0„cos0„ 


sin 02 
cos 02 
sin0„cos0„ 


2sin0^ cos 9^ 
-2sin0^ cos 9^ 
cos 02 -sin 02 


The Q matrix can then be calculated from equation (241): 


Q =T-^Q RT R-^ 

n n 


(240) 


(241) 
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The extensional matrix A can now be calculated from equation (242): 


plies 

^ Qkh = 


k=l 


l.lOxlO^ 

3.70x104 

0 


3.70x104 

l.lOxlO^ 

0 


0 

0 

3.65x104 


M 

in. 


( 242 ) 


The coupling stiffness matrix B, for a nonsymmetric matrix, can also be calculated from equation (243): 


plies 

^ ~ Qk^k^k 


k=l 


1.11x10^ -1.11x10^ 0 

-1.11x10^ 1.11x10^ 0 

0 0 -1.11x10^ 


M 

in. 


(243) 


The bending stiffness matrix D, equation (244), completes the calculations of the stiffness terms: 


plies 

D= 12, 

k=l 


^,3 ^ 

iL + t 
19 ^h^k 


5.73x102 

1.93x102 

0 


1.93x102 

5.73x102 

0 


0 

0 

1.9x102 


M 

in. 


(244) 


The global stiffness matrix from equation (246) can then be determined by combining the A, B, and D 
matrices: 


Global = 


l.lOxlO^ 

3.70x104 

0 

3.70x104 

l.lOxlO^ 

0 

0 

0 

3.65x104 

l.llxlO^ 

-l.llxl03 

0 

-1. 11x103 

l.llxl03 

0 

0 

0 

-l.llxl03 


l.llxl03 

-l.llxl03 

0 

-l.llxl03 

l.llxl03 

0 

0 

0 

-l.llxl03 

5.73 xl02 

1.93x102 

0 

1.93x102 

5.73x102 

0 

0 

0 

1.9x102 


.(245) 


Equation (246) is written in terms of the global stiffness matrix as 




^Ox 



^Oy 


= [Global] 

^Oxy 

Kx 

My 



"xy 


^xy 


( 246 ) 
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Additionally, equation (246) can be solved to determine the effective material properties 
of the composite beam: 


E^ = -^ = 3.41x10^ 


he. 


Ox 


m. 


V = = 0.42 


= 7^ = 3.41x10^ 

'■ " in.^ 


he. 


Oy 


V = = 0.42 

^ £0y 


s Ibf 

G^, = — ^ = 1.20x10^ 


hv — ,_^2 
'‘■I oxy 


in. 


The analysis was repeated for each of the three composite beams. Table 12 contains the effective 
material properties for each beam. 


Table 12. Effective material properties. 



8 

(Ibf/in.2) 

(lbf/in.2) 

^xy 

(lbf/in.2) 

Beam 1 

[90,90] 

4.12x105 

4.12x105 

7.50x10^ 

Beam 2 

[45,45] 

2.40x105 

2.40x105 

1.88x105 

Beams 

[90,45] 

3.20x105 

3.20x105 

1.15x105 


(247) 


6.5 Comparison of Test Data to Analytical Predictions 

The effective material properties can now be used with the elastic beam equation (124) to calcu- 
late the deflections in the beam. To illustrate the steps in the analysis, the results obtained for beam 1 
at a typical load increment follow: 


f2(^) 


J_ 

El 


{Pax^ + C^x -I- 



(248) 
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E = 4.12x105 

in.^ 

(249) 

bh^ (I.189in.)(0.265in.)^ 4 

= = 0.00184 in.4 

(250) 

12 12 


X - 4.375 in. 

(251) 

C 3 = -25.39 Ibf/in.^ 

(252) 

C 4 = 8.707 Ibf/in.^ 

(253) 


^2 


(1.94 lbf)(3 in.)(4.375 in.)" 


(4.12X105 ‘“1 

f0.00184 in.4] 

i Ibf 1 

1/ \ 1 

( Ibf ^ 

V in.^ } 

\ / 

+\ -25.39^ 

4.375 in.) + 

00 

0 

( 

0 ^ 



V inEJ 


V m.'’ J 


(254) 


J 2 = -0.062 in. (255) 

The above analysis is repeated at incremental loads to generate a load versus deflection curve for each 
beam. 


The results of the three-section bending test can now be compared to the results calculated using 
the elastic beam equation with the effective material properties derived from the rule of mixtures and the 
laminated composite plate theory. A measure of the possible error associated with these calculations is 
determined as follows: 

• Error from Young’s modulus test = 6.4 percent 

• Error from shear modulus test = 5.6 percent 

• Error % = ^|i6A%f + (5.6%f = 8.5% . 

Error bars on figures 82-84 will reflect this measure of error. 

As shown in figure 82, the predicted deflections for beam 1 [90,90] are in good agreement with 
the test results, indicating that the effective Young’s modulus used in the predictions, = 4.12 x 10^ Ibf/ 
in. 2 , adequately predicts the behavior of this two-layer composite. 
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Figure 82. Beam 1 comparison of predicted versus test deflections. 
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Figure 83. Beam 2 comparison of predicted versus test deflections. 
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Figure 84. Beam 3 comparison of predicted versus test deflections. 
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Figure 83 compares the predicted deflections for beam 2 [45,45] with test results. The data 
compare reasonably well with the test results, indicating that the effective Young’s modulus used in the 
predictions, = 2.31 x 10^ Ibf/in.^, also adequately predicts the behavior of this two-layer composite. 

Figure 84 illustrates the results for the nonsymmetric beam 3 [90,45]. The test data for beam 3 
[90,45] compare very well with the test results. The effective Young’s modulus calculated for this beam, 
= 3.2 X 10^ lbf/in.2, was calculated using the same methods as the previous symmetric beams. 

Based on these results, it can be concluded that the solution to the constitutive equations pre- 
sented can be used to accurately determine the effective material properties for a symmetric or a 
nonsymmetric multilayered laminated composite. 

6.6 Analysis of Beams Using Finite Element Methods 

The analysis of simple composite stmctures such as beams can be achieved using the methods 
described in this TM. However, when the geometry or loading of the structure becomes more complex, 
this method of analysis may be too difficult to apply. For that reason, it is important to extend the verifi- 
cation process to include finite element modeling. 

The input for the MSC/Nastran program is calculated from the equation for the stiffness matrix, 
equation (257), = [5„]~^. The MAT2 input for beam 1 can be taken from the matrix. For example, 

the stiffness matrix for layer one of beam 1 [90,90] is calculated and given: 




4.30 xl05 
6.45 xlO^ 
0 


6.45 xlO^ 
4.30 xlO^ 
0 


0 

0 

7.50 xl04 


(257) 


Incorporating these values into the MAT2 results in the following: 

$ Description of Matrix Material: 

MAT2, 1,4.30+5,6.45+4,0.0,4.30+5,0.0,7.50-1-4, 

0., 0., 0., 0., 0. 

The orientation of each layer is described in the PCOMP command. An example of the PCOMP 
command for beam 1 is provided below for reference: 


$ 

PCOMP, 1, , , 10000.0, STRN, 

+, 1, .1325, 90.0, YES, 1,.1325, 90.0, YES 

$ 


The finite element model is generated in the same procedures as those verified in section 4.12. 
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Figures 85-87 show that the results of the finite element models for each of the three beams 
are in agreement with both the test data and the analytical solution from the elastic beam equation (86). 


— Beam 1 [90,90] — * — Beam Equation — Q— pjpjte Element Model 



Figure 85. Beam 1 comparison of finite element versus test deflections. 


-♦ — Beam 2 [45,45] Beam Equation B Finite Element Mode 



Figure 86. Beam 2 comparison of finite element versus test deflections. 
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Figure 87. Beam 3 comparison of finite element versus test deflections. 


6.7 Summary and Conclusions 

The preceding results indicate that cementitious composites can be analyzed using the laminated 
composite plate theory. This is a significant improvement over the methods based on the rule of mixtures 
and the transform section methods. The composite plate theory has no limitations on the orientation of 
the fibers or the geometry of the layup. As a result, the designer can strategically place and orient the 
reinforcement fibers to optimize the design. 



0 0.02 0.04 0.06 0.08 0.10 0.12 

Deflection (in.) 
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7. CONCLUSIONS AND RECOMMENDATIONS 


7.1 Conclusions 

This research establishes the procedures for producing, instrumenting, and testing graphite- 
reinforced cementitious specimens to determine the material characteristics. Tensile tests are used to 
successfully determine Young’s modulus and Poisson’s ratio. losipescu shear tests are used to determine 
the shear modulus. Comparisons are then made to the material properties calculated from the mechanics 
of materials approach for composite analysis; i.e., rule of mixtures. Consistent with the literature,^’^’^^ 
it is determined that the rule of mixtures is accurate for determining the effective material property, E-^, 
provided the fibers are continuous and aligned in the direction of the load. However, the remaining 
material properties, £’2, v^2’ G^2 accurately predicted by the rule of mixtures. 

A modified transform section approach is used with the £j value from the rule of mixtures to 
accurately predict the deflection of a beam in pure bending with a single layer of graphite reinforcement. 
A similar approach is used for a multilayered beam in pure bending. The beam analysis is extended to 
include nonlinear material properties. The results compare well to test data. 

Laminated composite plate theory is investigated as a means for analyzing even more complex 
composites, consisting of multiple graphite layers oriented in different directions. Equations for the 
effective material properties are derived based on the laminated composite plate theory. The equations 
are arranged to form a system of simultaneous equations from which an exact solution can be found. A 
Mathcad solution sheet is developed to make the calculations. An increase in accuracy is demonstrated 
when this solution method is compared to other solution methods. The results of the Mathcad solution 
sheet are verified with an example problem. A composite plate with brown material properties is ana- 
lyzed using the Mathcad solution sheet. The plate is then analyzed using a finite element model. The 
results compare very well, demonstrating the accuracy of the solution method. 

The analytical methods developed in this TM are demonstrated with the analysis and test of three 
beams with various ply angles and layups. The effective material property for each beam is determined 
by using laminated plate equations and the Mathcad solution program. In addition, a finite element 
model is developed for each beam. The results of three-section bending tests are then compared to the 
analytical results to determine the accuracy of the methods developed in this TM. The analytical results 
compare very well to the test results. 

The research presented establishes a methodology that can be used to determine the deflections 
and stresses in graphite-reinforced laminated composites. In addition, it demonstrates methods of deter- 
mining the appropriate material properties necessary for the analysis. 
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APPENDIX A— MATHCAD SOLUTION SHEET: RULE OE MIXTURES 


Note: This program is provided as an example only. The content of this file is similar to that 
found in the text. However, this program is subject to modification and may not match the previous text. 

A.l Rule of Mixtures for a Composite Material 

The rule of mixtures is a method for determining the effective properties of a composite material 
based on the contributions of the individual constituents. This method assumes that the strains in the 
fiber direction of a composite are the same in the fibers as in the matrix. Since the strain in the fiber 
is the same as the strain in the matrix, the sections normal to the fiber direction remain plane after 
stressing: 

• Determination of 






At 






The strain in direction 1 , the fiber direction, is given by 

AL* 



The strain in the fiber direction is same for the fiber and the matrix: 

*^fiber . 

Applying Hooke’s law and assuming both the fiber and the matrix are linear: 

cfinatiiK . 
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This relationship is similar to two springs acting in parallel. The total load (P) is carried by the load 
in fiber (Py) and the load in the matrix (P„): 


P- Pf T Pm' . 


"Wl/VH 


-V144^ 



The force equation can then be written in terms of the stress and area: 

aj Ai := ofiber'AfJber + •^rtiatrix'Arriatrix • 

Now substitute the Hooke’s law equations for the stress terms to relate the force equation in terms 
of strain and Young’s modulus: 


Epsi Ai E^er-si-Afiber + Ertiatiix'^l ■^matrix • 


Simplifying the equation, 


Ptiber’^fiber Pniatrix' ^matrix 

El := + 


Al 


Al 


The fiber volume fraction and the matrix volume fraction (Yj^atrix) defined as 


^matrix •“ 


Amatiix 

Al 


Afiber* 
^fiber:* — . 

Al 


In addition, the fiber volume fractions must satisfy the equation: 


Afiber + ^luatnx ■- 1 • 
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Thus the expression for the effective Young’s modulus in direction 1 is as follows: 


El ^flber'Eiiber . 

Determine the effective property in direction 2, perpendicular to the fiber direction, 



Deformed 


The fiber and matrix transverse strains are not assumed to be equal. However, in the transverse 
direction, the cross-sectional areas are equal and the load is constant across the area; therefore, the stress 
is equal for the fiber and the matrix. Then the strains in the fiber and matrix can be written in terms 
of the stress: 


and 


®fiber •- 


■ 

P2 

Efiber 



I 

^matrix •- Z . 
Ematiix 


The total transverse deflection is the sum of the fiber deflection and matrix deflection: 

Aw ;= Awj5ber-l- 
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Defining the strains in terms of the deflection: 


Aw* 



imer* 

^fiber “T , 
Miber 

and 

Awn^trix 

^matrix . 

Mtiatrix 


The transverse dimension over which the fiber and matrix strains act can be written in terms 
of the fiber volume: 


and 


Lfiber := ^fiber w"* 
hmatiTX •= ^matrix' ^ • 


The deflection equation can now be defined in terms of strain: 

S2'W':= sfiber(Vfiber-w) + SKatiix(VmatrixW')*- 

Now the w term vanishes from both sides of the equation. Hooke’s law equations 
can be substituted for the strains: 


02 02 02 * 

~ •= '^fiber r + ^matrix’ r . 

^2 ^hnatiix 


The stress term now vanishes from both sides: 


1 _ VfJber ^ ^matrix 

^2 Efjjer Ematiix 


Solving the above equation for the E2 term. 


EfiberTmatrix 

E2 := 

^matrix' Efiber + ^fiberTniatrix 
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Determination of the Poisson’s ratio between the transverse and the fiber directions: 
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Determine the shear modulus (G12) by assuming the shearing stress on the fiber and the matrix 
is the same. This assumption means that the shear deformation can be written as 

I 

T 

^matrix T . 

’^matrix 

• Shear strain in the matrix: 


• Shear strain in the fiber: 


tfiber := 


T 

Gfiber ■ 


712 := 




Shear strain in the fiber and transverse direction. 

The shear deflection is a function of the shear strain and the width (w): 


and 


A := 7' w 


A 


fiber + ^matrix: • 
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As shown in the previous equations, the deflection for the matrix and fiber can be stated in terms 
of the volume fractions, and V^i^an-ix- 


and 


^matrix •“ fnriatix' (^matrix' 
ifiber := 7fIber(Vfiber- w)*. 


Then, substituting into the deflection equation and allowing the w term to vanish: 

y Tmatrx' ^matrix + ffiber'l^flber • 

Finally, substitute the Hooke’s law equations for shear strain and allow the shear stress: 


— := Vmatiix T + ^fiber T . 

*^12 *^Kiatrix Gfiber 


Now, this equation can be simplified by allowing the shear stress term to vanish and solving for G^2- 


Gi2 


*^Hiatrix*^fiber 

Vmatrix'Gfiber + l^fiber Gmatrix ’ 


Calculate the effective material properties of a cementitious composite: 
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Material properties: 


Ematiix ■- H5-10 


Ef5ber:= 33.510" 


Vnifltriy ■— 0.27 


vfiber:=0-30 


Pinatm ■“ 35- 


P fiber - 064 - 




^matrix 

3(1 + ^matrix) 


Gfiber - 


^matrix ~ 4.528 x 10 


Gfiber= 1-2S8X 10 


Fiber area in a graphite yam: 


Ayani — •000175-in 


Nyam — ® 


Afiber ■“ Nyam'^yam 


Aj := .12 in^ 


^flber ■- 


^matrix - 1 “ ^fibei 


Vfiber = 0-012 


^matrix = 0.9S8 . 


Thus, the expressions for the effective material properties can be calculated: 


Young’s modulus in direction 1: 


E^l •- ^fiberTfiber + ^rnattix Ematrix 


El = 5.045 X 10^ — . 

.2 

in 
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Young’s modulus in direction 2; 


and 


^ 2 -= 


Efiber'Enriatrix 


^matrix' Efiber + ^fiber'Ematiix 


E2 = 1.164 X 10 


5lbf 


• Poisson’s ratio: 


and 


^12 := ^ fiber Yfiber + '’matiixYmatris 


V12 = 0.27 . 


• Shear modulus: 


G12:= 


•^matrix' ^fiber 

Vmatrix'Jfiber + Vfiber'Gmatrix 


and 


4lbf 

Gi 2 = 4.581x10 — . 

in 


• Composite Density 


Pcomposite ■“ PfiberYfiber+ P matrix ^matrix 


and 


Pcomposite ~ 34.592 ■ 


Ibf 


Now compare the results from the rule of mixtures to the results of tensile testing. 
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^0,1 

Strain 


Side 1 
0-© Side 2 
Average 


Young’s modulus from the tensile test: 


(1500 - 500) 


Ibf 


Elest •= 


(3000 - 100Cf)10 


5 Ibf 

^test = 5 X 10 - 

in 


A error := 


Etest-El 


Etest 


■10 


id 


Aerror = 0.898 • 



Krenchel methods: 


• Calculate the effective material properties of a cementitious composite: 


Efiber ■“ 33.5 10 ■ ^ 

in 


3 Ibf 

Ematiix ■“ 115- 10 ■ - 
in 


:= 0.27 vf5b« := 0.30 


... Ibf 

Ibf 

P fiber - 064- — 

Pmatrix — 35- 



Efjjg]. 





Gfibei:= 1.28846 X 10^ — 
in 

• Fiber area in a graphite yam: 

2 

Ayani := .000175-in 

^flber NyanvAyani 
wbeam := 1-1^3-in 

Ai = 1.40141 X 10“ 


*^matrix “ 4.52756 x 10 — - 

A 

m 


Nyarti ® 

0.241 -m 
tpljr^=^— 

^1 •= %ily'^eam 


Afiber 
Vfiber := — — 
Al 


^matrix 1 ~ ^fibei Vfiber = 9-9899 x 10 


^matrix ~ 9.9001 x 10 


Thus, the expressions for the effective material properties can be calculated: 

El ;= ^fiber'EfJber + ^jmtrix'Ematrix- 
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Young’s modulus in direction 1: 


El =4.48313 x ; 10^ M . 

.2 

m 


• Young’s modulus in direction 2: 


Efiber'Eifiatiix 

E-2 := 

^matrix' Efiber + ^fiber ^matrix 

and 

£2=1.16156x10^— . 

. 2 
m 


• Poisson’s ratio: 


and 


^12 := '' fiber Yfiber + '’matrixYmatrix 


V12 = 2.703 X 10 


• Shear modulus: 


and 


• Composite density: 


G12:= 


*^rtiatiix*^fiber 

Vnratrix'^fiber + ^fiberG’matrix 


Gn = 4.57308 x 10^— . 

. 2 
m 


and 


P composite ■“ P fiber YfJber + P matrix Yniatiix 


P composite - 3.4651 x 10 


1 M 

. 3 
111 


Now refer to Krenchel, p. 29: 


r|/:=45 
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• Krenchel case 2a 

ai := 1 


• Krenchel case la 


F actor 



Factor = 3.5 x 10 


Gfc (^rnatiix '^Kiatiix) 


+ 


^fiber'Efibe 



Gfc= 1.28489x10^— . 

. 2 
m 


cos u/ 

V uso ^ 


2V 

A 


s-2 := 1 


n := 1 P := Vfjbgi 


k:= 


1 


1-^ 

2 

''matrix 


k= 1.06786 X 10 


''matrix 

vcomposita := Vcomposita = 2.67498 ^ lO' 

M2 := [k + ai n p) - (k- Vi[,attix '’corftposite) W2 = 1. 00072 x 10^ 

Efc •= ®2 Emattix 


Efc 


= 1.15083 x 10^— . 

. 2 
m 


Efc ■- _Ejfia{ji)t k (^l - + Ef5j3gj 


Eft = 4.48513 X 10 


5M 
. 2 


m 
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42:^ 













► 




1 






Direction of the Load 


Krenchel case 2b 


:= 0 

180 


uro := — --^O 
180 


These angles describe the load in relation to the fibers (see KrencheP, pp. 17 and 26). 
The n represents the number of groups of fibers; in this case, two groups: 


n := 2 ♦ 


This represents the proportion fibers in each group; half of the fibers are in group 1 and the other half 
are in group 2: 


Efficiency factor: 


ai := 0.5 a2 := 0.5 


r| := ai cos[i|/i] + a2 cos[i|/2) 


ti=5 3<10 ^ • 


Makes the derivations easier (p. 27): 


k:= 


1-P 

, 2 
i '’matrix 
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Page 34 equations: 


Volume fraction for the fiber: 


'’matrix + (ti-.i-np) 


'’composite ~ 2.73399 x 10 


P := Vfjbei 


f, = 9.9899 X 10 


Guess values to make the program solve: 


5x:= 110 


-1 


5y:= 110 


5 Ibf 

^composite 1010 ? - 

in 


Equations (28) and (32) from p. 27: 

Ematrix'k (sx + '’matrix'^y) 

^composite - 


Efjber p cos(\|/i) •[(5x) cos{\(;i) + Sysin{r|/i)^ 
\ / \2 . / 


1 





Equations (29) and (33) from p. 27: 

0.0 = [En^trixk'(''rfiatrix'®x + ^y)] + 


2 2 2 
Efiberli sin('Fl) + 5ysin('|Jl) _ ■ 

2 " 2 2 
+ EfiberPcos(i|/2)^[(5x)cos(i|/2) + Sjr sin(i|/2)'' 


Equation (34) from p. 27: 


-By = Vccjjfippsite'Ex • 



Solve three equations with three unknowns: 


^composite '' 

By •= Find[Ecomposite>®y.5x) 


5 Ibf 

^composite ~ 4.48513 x 10 ^ 

in 


Now use equation (39) to calculate shear modulus, p. 29: 


■^composite 


:= (l - P) GmatrK+ _P Efiber (ai sin(i|/l)^ cos(i|;i)^ + a2 sm(i(;2)^ cos(i|/2)^)_ 


'^composite ~ 4.48233 x 10 
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Axial Stress (psi) Axial Stress (psi) 






Microstrain (;u in./in.) 


Figure B.4. Average transverse strain Dec_28_4. 










Axial Stress (psi) Axial Stress (psi) 
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Axial Stress (psi) Axial Stress (psi) 


y=0.398x + 16.653 



Figure B.9. Average tensile strain Dec_28_7. 


y=-3.0643x +3.6812 



Microstrain (/^in./in.) 

Figure B.IO. Average transverse strain Dec_28_7. 
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Shear Stress (psi) Axial Stress (psi) 



Figure B.13. Average tensile strain article SG4. 


Shear 4 Shear 2 



Figure B.14. losipescu shear test articles. 
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APPENDIX C— MATHCAD PROGRAM FOR COMPOSITE BEAM IN BENDING 


Note: This program is provided as an example only. The content of this file is similar to that 
found in the text. However, this program is subject to modification and may not match the previous text. 

C.l Mathcad Solution Sheet: Calculate Deflection for Three-Section Bending Test — One 
Graphite-Reinforced Section on Top With Linear Material Properties 

C.1.1 Calculate the Deflection for a Beam in Pure Bending 

The rule of mixtures method for calculating the effective material properties for a composite 
beam will be verified by comparison to test data. A composite beam was made to the following 
dimensions: 

• Graphite strands: 


Wg:=1.20in hg:=.115in . 


• Composite section: 


he := .361in- hg We:=1.20 in he = 2.46 x 10 ^ in • 


I Graphite Section 


Matrix Section 


J 



The fiber volume can be determined from the geometry of the test specimen. This test specimen 
has nine yams across the width: 

2 

Ayam ■= 000175-in 
Aj := .115-in- 1.20-in 


Nyam ■“ ^ 


^fiber: 


Aiiber 

Al 


Afiher-- Nyam'Ayajti 
Vf5ber= 1-141x10'^ 
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The volume fraction for the matrix: 


^matrix ■- ^ “ ^fibei ^matrix = 9.8S6 x 10 


The material properties of the graphite and matrix are the same as previously defined: 


3 Ibf 

^matrix ■“ 11510 - 

in 


Efiber ■“ 33-510 ■ ^ 

in 


^matrix ■- 0.27 


Vf5ber := 0,30 


P matrix ■“ 35- 


P fiber - 0b4 - 


'^matrix •- 


^matrix 

2'(l + ^matrix) 


Gfiber - 


2-[l + vfiber) 


4 Ibf 

‘^itiatiix “ 4.528 x 10 - 

in 


Gf5ber= l-288 x 10^ — 
in 


The effective material properties for the composite section of the beam can be calculated: 


Young’s modulus in direction 1: 


Young’s modulus in direction 2: 


Poisson’s ratio: 


^1 •” ^fiber'Efiber + ^matrixTmatrix 


El = 4.9602 X 10 


EfiberTmatiix 

^matiixTfJber + ^fiberTmatrix 


Et = 1.163 X 10^ — 
. 2 
m 


V12 - ''flberYfiber + '’matiixTWtrix 


V12 = 2.703 X 10 
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Shear modulus: 


G12:= 


^matrix' *^fiber 


^matm-Gflber + ii^f3bert>matrb: 


Gi 2 = 4JS»10 


4lbf 


• Composite density: 


P composite •“ P fiber' ^fiber + P matrix' ^matrix 


P composite 


3.46 X 10^ — 
. 3 
in 


Now consider the composite beam composed of two homogenous materials. The top section of 
the beam is made from a material whose properties are calculated from the mle of mixtures. The remain- 
der of the beam is made from the matrix material. The moment of inertia for the beam can now be 
calculated using the transform section method: 

El 

E^trix ‘ 


• Ratio for transform section: 


p = 4.313 . 

• Transform section by multiplying the width of the graphite sections by rj: 

Wg := r| Wg 


Wg = 5.17(5 in . 

•Now calculate the centroid of the transformed section: 


Ag:=hg Wg . 


• Calculate the area of the top section: 


yg := he + 


th] 
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Locate the centroid of the top section: 


• Calculate the area of the matrix: 


Ac := hc'Wc. 





• Calculate the centroid of the matrix: 


Ybar:= 


Ag yg + Ac-yc 

Ag + Ac ■ 


Equation to find the centroid of the transform section: 


Ybar = 2.4366 x 10 ^ in 


C.1.2 The Centroid of the Transform Section 

The cross section of the composite is symmetric in the J direction. The x location of the centroid 
is the center of the beam: 


Wc 

Xbai-:-— Xbar=0.6in . 

• Calculate the moment of inertia for the top section of the beam by using the parallel axis theorem: 

b := Wg and h := hg . 

• Moment of inertia of the top section in respect to its centroid: 


• Area of the top section: 


Ig := — b 

^ 12: 


Ag:=b b . 
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• Distance from the center top section to the centroid of the section: 

d:=Ybar- 

• Applying the parallel axis theorem: 

Ig_total ~ Ig Ag- d • 

• Total moment of inertia for the top section with respect to the centroid of the entire beam: 

-3 4 

^g_total = 2.7874 x 10 in • 

Now, using the same methods, calculate the moment of inertia for the matrix (bottom) section 
of the beam: 

b := Wc 

and 

h:=hc, • 

• Area of the matrix: 

Ac:=b h . 

• Distance from the center of the matrix section to the centroid of the section: 

fhc') 

(1-Ybar- — . 

V 2 y 

• Moment of inertia for the matrix with respect to the centroid of the beam: 

2 

total “ + Ap d 

-3 4 

Ip total = 5.7864 X 10 in • 

• Total moment of inertia for the transform section: 

hotal Ic^total + Ig_tota] 

Itotal = S.5739 x lO^^in*^ • 
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The beam was then subjected to a three-section bending test. A rigorous analysis of the beam test 
was performed to establish a base of comparison to finite element modeling and test data. The compari- 
son will provide insight into the accuracy of the methods used for the analysis. 

• Free-body diagram of the beam: 



P 

◄ a ► 

\ 

M b ► 

\ 

P 

M a ► 







P 


P 


Derive the equation for the elastic curve from sum of the moment at x: 
• Section \ ,x<a: 


CT a 

El-^yi 

dx 


iPx 


ElAyi = --x^ + Ci 
dx 2 


P 3 

Elyi = - X +C1X+C2 
6 
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Section 2,x> a but x< b: 
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Summarize the equations. Replace the derivative with ft 


P 2 

EI01 = --^ + Cl 


Ei e2 = P ax+C3 


P 3 

Elyi = — X +C 1 X+C 2 
6 


P 2 

E'Iy2=— ax +C3X+C4 


P 2 * 

E 103 = 2 P a x + P b x - - X + C 5 


2 P 2 P 3 

EIy 3 = P ax H bx x +C 5 X+C 5 

2 6 


Values from the test fixture and the test loads: 


Maximum load from 5/2/00 test: 


Load =50 P := .2207-lbf- 


Load 


Geometry of the test fixture: 


a := 3.0 in b s? 2.75 in 


L:=8J5in 
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• Input values to represent the boundary conditions: 

XI := a X2 := (a + b) X3 := L 

Xq := O.O in xj = 3 X 10 ^ in X2 = 5.75 x 10 *^ in X3 = S .75 x lo'^in 

• Guess values to start the program to solve for the constants: 

Cl l.O lbf in^ C2 := l.O lbf in^ C3 := l.O lbf in^ 


C4 := l.O lbf in^ C5 := l.O lbf in^ := l.O lbf m^ • 

The elastic curve equations written to satisfy the boundary conditions: 

• For the boundary condition x = 0, at x = 0, yj =0: 

O.O lbf in^ = - Xo^ + Ci-Xo + C2 • 

6 

• For the boundary condition x = a, at x = a yi =y 2 '. 


P 3 P 2 

— XI + Ci'Xi + C2 = — axi + C3X1 + C4 
6 2 


P 2 

— XI + Cl = P-axi + C3 


For the boundary condition x = (a + b), at x = a + b y 2 = y^'- 


P 2^ 2P^2P3„ 

— ax2 + C3 X2 + C4 = P ax2 H b-x2 X2 +05 x2 + 0^ 

2 2 “ 6 “ 


P 2 

P a x2 + C3 = 2 P a x2 + P b x2 - “ + *^5 


For the boundary condition x = L, at x = L, y^, = 0: 


O.O lbf in^ = P a-3^^ + — b-:g^ - — X3^ + C5 X3 + Cg 
2 6 
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Find constants to solve the above equations simultaneously: 


C2 

C3 

C4 

C5 


Find(Ci,C2,C3,C4,C5,C6) . 


Ends the program loop. 

The above equations are solved simultaneously to find the constants. The constants are listed 

below: 


Cl = -7,681 X 10^ Ibf ir? 
C5 = -2.641 X lO^lbfin^ 


C2 = 0 X 10°lbfin^ 
Cg = 3.222 X lO^lbfin^ 


C3 = -1.169 X lO^lbfin^ 
C4= 4.007 X 10^ Ibfin^ • 


The value for the moment of inertia is set to the value calculated from the rule of mixtures 
and the transform section method. Young’s modulus for the matrix is used for the beam equations: 


1 •“ hotal 

1 = 8.574 x10 


E= 1.15 X 10^ ~ 
2 


in 


Restate each elastic curve equation as a function of x: 
• Section 1: x < a: 


yi(x) := ^ + Ci-x+ 02^ 
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Section 2: x> a but < a+ b: 


J_/L 

EI) 1,2 


y2(x) := - ax + C3 X+ C4 


Section 3, x> a + b: 


®2C=t) := :^(P a-3c +C3) • 
Ji-1 


If P 2 

03(x) := 2 P a x + P b x x + C5 

Ell 2 ■ 


if 2 P 2 P ? 

V3(x) := P a x H b x x + C5 X+ Cg I 

Ell. 2 6 ' • 


Check boundary conditions; x = 0 and yj = 0 and yj = y 2 at a" = a: 


yi(O.O-in) = 0 X 10 ^ in 


yi(a) = -1.93 x 10 ^ in 


0lCa) = -3,725 X 10 


-2 


y 2 (a) = -1.93 x 10 ^ in 


e 2 (a) = -3.725 x 10 


-2 


• lfx = a + b, then y 2 = y 3 ; 

y2(a + b) = -1.93 X 10 ^ in 02(a + b) = 3.725 x 10 ^ 

y3(a + b) = -1.93 X 10 ^in 03(a + b) = 3.725 x 10 ^ 

• Ifac = L, theny3 = 0: 

y3(L) = 0 X lO^^in 

C.2 Mathcad Solution Sheet: Calculate Deflection for Three-Section Bending Test — Three 
Graphite-Reinforced Sections With Nonlinear Material Properties 

C.2.1 Determination of Transform Section Properties 

Cross section of composite beam: 

• Number of layers: 


150 




• Depth of each layer (bottom surface is zero): 


D := 


^ .025 ^ 
.0875 
.150 
.2125 


\215Q) 


• Modulus for the graphite: 


E2 := 5.04- 10 


J 


• Modulus for the matrix 


El := 115 10^ • 


• Ratio for transform section: 


f] ■= 


El 


^ = 4.383 • 
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Dimension of the composite section 


he -= .075 ; wc := 5.00 . 


Dimension of the graphite section: 


Calculate the area of the fiber: 


Calculate the area of the matrix: 


lig := .050 Wg := 5.00 


Wg := r| 


Afiber •“ hg (wgj 


^matrix •- hc'^'c 


Equation to find the centroid of the transform section: 




Aman D; 




^bai:= 


i= 1 




'map 


i= 1 


^ Afiber 
^matrix 
Afiber 
^matrix 
\ Af 5 ber / 


Centroid of the transform section: 


Ybar= 1 . 5^10 


Moment of inertia of the transform section using the parallel axis theorem: 


I:=Ig+Ag-d 


T - ^ 


T ^ V, 3 
in := — Wr he 
c 12 " ^ 
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• Determine the distance from the center of the graphite to the centroid of the section: 

Layer_d(n/) for iel. ij; 

d- ^ Ybar- D. 

d 

d := Layer_d(\i/) 

^ 0.125 
0.063 
d= 0 

-0 063 
^-0.125 

• Moment of inertia for each section: 

I:=Ig+Ag.d^' 

Layer_I := for ie 1,3,5 

Ii<e-Ig+ A]^p. (d^^ 
for j € 2,4 

2 
I 

Ilayer := LayerJ 

0.017 
0.002 
0 

0.002 
0.017 

• Moment of inertia of the five-layer composite: 
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C.2.2 Calculate the Deflection for a Beam in Pure Bending — Consider the Nonlinear Nature 
of the Material Properties 

Beam is in pure bending (no shear) in section b. Free-body diagram of the beam: 



• Sample values from the test fixture and the test loads: 

t := 0.30 in Thickness of the beam 
Yloadj^lbf 

P := Load applied to beam . 

• Geometry of the test fixture: 

a:=5.25 in b := 4.5 in L :=15.0 in • 

• Input values to represent the boundary conditions: 

XI := a X 2 := (a + b) X 3 := L 

Xq := O.O in xj = 5.25 X 10*^in X 2 = 9.75 x lO^iti X 3 = 1.5x10^ in • 

• Guess values to start the program to solve for the constants: 

Cl := l.O lbf in^ C 2 := l;Q-lbf-in^ C 3 := l.O lbf in^ 

04 := 1.01bf in^ CyteWlbf-in^ Cg := l.O lbf in^ • 
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Following are the elastic curve equations written to satisfy the boundary conditions: 


For the boundary condition x = 0 and = 0, 


O.Olbfin^ = - Xo^ + Ci-Xo + ^2 
6 


For the boundary condition x = a; and = 62'. 


P 3 P 2 

— XI + Cyxi + C2 = — axi + C3X1 + 
6 2 


P 2 

— XI + Cl = Paxi + C3 


For the boundary condition x = (a + b); )>2 = y3 and 6^ = Oy. 


P 2 2 P 2 P 3 

~ a x2 + C3 X 2 + Cq = P a x2 H b x2 X2 + C5 X 2 + Cg 

2 2 6 


P 2 

P a x2 + C3 = 2 P a x2 + P b x2 - — X2 + C5 


For the boundary condition v = L; y3 = 0: 


3 2 P 2 P 3 

O.O lbf in. =P ax3 H b x3 X3 +05 x3 + 0^ 

2 6 


Find constants to solve the above equations simultaneously: 


C2 

C3 

C4 

C5 

yC6^ 


: Fmd(Ci,C2,C3,C4,C5,C6) 


Ends the program loop. 
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Now compare these equations with actual test results: 

• Values used in the beam equations: 

E= tl5x 10^— _7 4 

■ 2 I = 3.8 X 10 in • 

in 

The load is applied by hanging weight centered in the ‘b’ section: 

Test_load := P-2 • 

Deflection data at the center of the 15-in. beam {x = 7.5 in.) : 

y2(7.5 in) = -8.01449 x 10“ ^ in 
Testjoad = 5.67 x 10^ Ibf • 

Now repeat above analysis at each load step. Tload is the load data supplied to the elastic beam equation: 


^ 0 
6.5 
9.1 
11.8 
14.4 

17.9 

Yload = 25.0 
32.0 

40.9 
48.8 
56.7 
64.6 

J5.1 
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The deflection calculated by the elastic beam equation is contained in 


ycEii 


{ ° ) 

.092 

.129 

.167 

.203 

.254 

.354 

.452 

.578 

.690 

.801 

.913 


Vi. 06 ; 


This Excel file contains the test data for the five-layer composite beam in bending (the UAH 
canoe team supplied the test data): 


k:=2,3..89 


Test Data := 


Worksheet 


The test data is compared to the results of the analysis. 


Canoe Plate 5:^/.30 thick Super-7 



Test Data^ 4 

Elteflecticin in Center of Plate 

Test Results; 

Theory 
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This Excel file contains the results of tensile test specimens. These data will allow the program 
to consider the nonlinear properties of matrix material (tension). 



Worksheet 


m:= 2,3.. 69 


250 
200 
150 

w TestOjj^ 3 

100 
5i 
0 

• Assign values to stress and strain data matrix: 

*^concrete^ 3 ^conciete^ •= g • 

• Guess values to initialize the program: 

bo := 1 bi := 1 b2 := 1 n := 69 . 


Tensile Test of MatnxMatenal 
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Use Balaguru’s methods^ (least-squares methods): 


Given 


.2 ^ 


boin-bj-^ Pconeuete^ + ^2'^_ [<^conc]steJ ^concrete. 

i=l i=l i=l 


i' 

\2 ^ r ■|3 


bo- 5 ; *^coiicrete^ + I'^concreteJ + ^2'^^ |*^concreteJ - *^concrete.- ^concrete. 
i=l i=l i=l i=l 


2 3 4 " 2 

bo- ( °'corucreteJ ['^concrete J +^2'^ [*^coricreteJ [(^concrete J -^conciete 


i= 1 


i = l 


i=l 


1 = 1 


Find the minimum solution to the simultaneous equations: 

bi := Minerr(bo,bi,b2) 

Vb2^ 


The equation for strain from the resulting values: 


«p :- bo + bi'CTcoj^j^te + h2-Oconcrete 


Compare the results of the derivation with the actual test data. 



Cufve fit 

Test Data 




Now, increase the stress to levels exceeding the tension test: 


• Stress equation: 


= 5:87316 x 10 


2 Ibf 
. 2 


m 


Assume the stress varies from 225 to 800 in increments of 25: 


^ := 223,250.. 800 • 


• Equation for the ‘simulated’ stress values: 


s(^) :=bo + bi^ + b2^^ • 

Plot the simulated stress and the results from the tensile test. The simulated stress projects 
the curve to higher strain values: 



<p.E(y 


Curve fit 
Projected Data 


Young’s modulus for the matrix changes from the stress increases to 800: 


E= 1.15 X 10^ ~ 
. 2 
in 


^ := 800 




3(<f))l0 


-6 El = 8.07035 x 10 
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In section 2 for the elastic beam equation, x> a but <a + b: 


Ttloadj^ = 56.7; s:= 7.5 iD . 


The beam equation has been modified to calculate the deflection at stress levels above the linear 
range of the material: 


y2(x,4>) := 


f 


if 


<*> > 250,- 


V 


!( 4 >) 10 ' 


■,115000 

/ 


Ibf 


. 2 
in ) 


P 2 

, , — -a-x -H Ct-x + Ca 
^ 1 2 


72^7.5 m,4.23799 x 10 


-0.90774m • 


The deflection is calculated at increasing load levels: 


( ° 1 

0.092 

0.129 

0.167 

0.203 


0.254 


ycal_2 := 


0.3547 . 

.479 
.655 


.825 

1.009 

1.21 


V 1.49 ; 


Compare the results of the analysis with the results of the canoe panel-bending test: 
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GTano e Plate;; K bh^tine ar Prop eilie s 



Deflection in Center of Plate 

Test Results 
Non-Linear Theoiy 


The analysis is repeated for the properties from compression testing of the cementitious material. 


This Excel file contains data from compression cylinder testing of the cementitious material: 

Test_Comp := 

Worksheet 


i:=2,3..50 

• Equation to calculate the compressive Young’s modulus: 


Ecomp ■“ 


Test_Comp^j g 


Test^Compj^j ^10 


• The highest value is at kj = 44: 

Ecomp = S.12316 X 10^^ • 

The beam equation has been modified to calculate the deflection at stress levels above the linear 
range of the material: 

x:= 7 ;y in 


y2(^>EcQinp) 


Ecomp" 


Tbf 


P 2 
— ax + Ctx + C4 
.2 
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• Sample calculation at maximum load value: 

Y-jj 29616 x 10^) = -1.26312in • 

The results of the analysis are listed in_ygjjj_3: 


^ 0 
0.092 
0.129 
0.167 
0.203 


ycEilJ := 


.264 • 

.389 


.524 


.731 

.976 


Vi. 26 ; 


The results using the compressive data compare very well to both the canoe panel data 
and the results from using the projected tensile. 



Deffection in Center of Plate 

Test Results 

O O Non-Linear Tensionljata 

N on-Line ar Compressign D ata 
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APPENDIX D— MATHCAD SOLUTION SHEET: LAMINATED COMPOSITE 
PLATES— SYMMETRIC AND NONSYMMETRIC SOLUTIONS 


Note: This program is provided as an example only. The content of this file is similar to that 
found in the text. However, this program is subject to modification and may not match the previous text. 

D.l Calculation of Effective Material Properties 

It should be noted that for laminates that are symmetric about the midplane; i.e., the layers above 
the midplane are identical to the layers below the midplane, the B matrix will be zero. In the case of 
symmetric laminates, the equations for the effective properties can be greatly reduced. However, for the 
purposes of this analysis, the case of the nonsymmetric laminates will be considered: 




''All Ai2 Ai 3 Bii Bi2 B^^ 



Ny 


A21 A22 A23 B21 B22 B23 


SOY 

Nxy 


A31 A32 A33 B31 B32 B33 


TOxy 

Mx 


Bll Bi 2 Bi 3 Dll D12 ^13 


Kx 

My 


B2I B22 B23 D21 D22 D23 


Ky 

^Mxy j 


, B31 B32 B33 D31 D32 D33 ^ 


l^xyy 


This relationship is a system of six equations and six unknowns. The midplane strains and curvatures 
are the unknowns. The A, B, and D terms are derived previously and are now constants. The force 
and moment matrix can be defined to determine the effective material properties for the composite: 




^Aii-sox+ A12-50Y+ Ai370xy + Bu K^ + Bi2-Ky -1- Bi3-Kxy V 

Ny 


A21-50X + A22-50Y + A2370xy + B2lKx + B22% + B23'Kxy 

Nxy 


A 31 EOX + A32-50Y + A33-y0xy + B 31 KX + 832 % + B33-Kxy 

Mx 


Bii-eox+ B 12 BOY+ Bi370xy+ Dll-Kx+ D12% + Di3-Kxy 

My 


B 21 EOX + B22-S0Y + B2370xy + D 2 I + D22 % + D23-Kxy 

l^Mxy j 


i,B31-eox+ B32-50Y+ B3370xy+ D31-Kx+ D32 %+ D33-Kxy ; 


The solution to this system of equations can best be demonstrated through an example problem. 
The calculations for this problem will be made in a Mathcad Pro 2000 program file. Consider a lami- 
nated plate made of four layers of fiber-reinforced material. The angle of the reinforcements is 0° and 
45° and each layer is 0.005 in. thick. 
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Material constants are obtained by tensile and shear testing of the material: 


El := 20.01 10 


(5 M 

. 2 
m 


Young’s modulus in the principal material direction 1: 


E2:= 1. 301 10^ — 
. 2 
m 


Young’s modulus in the principal material direction 2: 


Gi2 := 1.001 10 


,6 Ibf 

. 2 
an 


Shear modulus in the principal material direction 12: 


V12 := 0.30 


Poisson’s ratio: 


V21 


V12-E2 

El 
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These properties are used to determine the compliance matrix: 
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The center of each ply must be located: 


*2 


L 4- 
1 2 




L V 


t- H 

3 2 


The transformation matrix is derived from the stress transformation equations: 


T := 
n 


cos 


(en)^ 


sm(eii)" 


sin (en)^ 2-sin (en) •cos (en) 

2 

(en)^ -2-sin(0n)cos (en) 


cos 


2 2 

-sin (en) ■cos (en) sin (er^ •cos (en) cos (en)^ - sin (en)^ 




^1 = 


1 X 10 Ox 10 Ox 10 


Ox 10 1 X 10 Ox 10 


, 0 ^ 


„0 


0 


„0 


vox 10 Ox 10 lx 10 ; 


^2 = 


■^3 = 


- 1 


- 1 

0 

5 X 10 

5 X 

10 

1 X 10 

- 1 


- 1 

0 

5 X 10 

5x 

10 

-1 X 10 

- 1 


- 1 

0 

-5 X 10 

5x 

10 

0 X 10 ^ 

- 1 


- 1 

0 'v 

5x 10 

5:X 

10 

lx 10 

- 1 


- 1 

0 

5x 10 

5 X 

10 

-1 X 10 

- 1 


- 1 

„0 

-5 x lO 

3 X 

10 

0X10 ^ 
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Reuter’s matrix: 


/ 0 0 0 '• 
1x10 OX: 10 0X10 

ox; 10° 1 ;^^ o:x;io° 

V O X 10° 0 X 10° 1 X 10°,; 


R 


A 0 0 ^ 
0 1 0 
1^0 0 2J 


The Qbar matrix represents the stiffness and transform matrix: 
Each ply has a Qbar matrix: 


^2.01 X 10^ 3.93 X 10^ 0 ^ 


Qbar^ = 


3.93x10^ 1.31x10^ 


. 2 
m 


Qbar^ = 


Qbar^ = 


Qbar^ = 


1 0 

0 

1 X 10 

J 


/ 6 

6 


5 ) 


6.56 X 10 

4.55 X 10 

4.7 X 10 



6 

6 

6 

Ibf 

4.55 X 10 

6.56 X 10 

4.7 X 10 


2 

6 

6 


6 

in 

[ 4.7 X 10 

4.7 X 10 

5.16 X 10 J 


/ 6 

6 

. 

6 'i 


6.56 X 10 

4.55 X 10 

4.7 X 10 



6 

6 

6 

Ibf 

4.55 X 10 

6.56 X 10 

4.7 X 10 


2 

6 

6 


6 

in 

4.7 X 10 

4.7;X 10 

5.16 X 10 > 


/ 7 

s 

o' 



2.01 X 10- 

3.93 X 10' 

0 X 10 



5 

6 

0 

Ibf 

3.93 X 10 ■ 

1.31 X 10 

0 X 10 


1 

0 

0 

6 

in 

ox io 

0 X 10 

1 X 10 , 

/ 
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Calculate the total thickness of the layup: 


plys 

h:= y t _2 

^ " h = 2 X 10 in 

n= I 


The A matrix or the extensional stiffness matrix is calculated: 

plys 
A:= ^ 

k= 1 


A = 


/ 5 4 4 \ 

2.67 X 10 4.95 X 10 4.7 x 10 

4.95 X 10^ 7.87 x 10^ 4.7 x 10^ 

V 4.7 X 10^ 4.7 X 10^ 6.16 x 10^ ^ 


Ibf 


The B matrix or the coupling stiffness matrix is now calculated. In the case of a symmetric 
laminate, the B matrix will be zero: 


plys 



E 

Qbair^ 

,W"k 

k=l 



^ 0 

-0 o '' 


-0 

-0 0 

Ibf 

kO 

0 0^ 



The D matrix or the bending stiffness matrix is calculated: 


plys 
k= 1 




D = 


( 1 
1.23 X 10 

6.09 X 10“ 

^3.92 X 10“ 


- 1 - n 

6.09 X 10 3.92 X 10 

131x10° 3.92 X 10“^ 
3.92x10“^ 1.01xl0°,j 


Ibfin • 
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To find a solution to this system of equations, define a matrix V to represent the forces and 
moments at the midplane of the plate. Then construct a “Global” matrix that represents the stiffness 
matrix. A computer subroutine “Isolve” will be used to find the matrix for the strain and curvatures. 
For the following equations to work in the Mathcad program, the matrices A, B, and D, must be made 
dimensionless: 


A:= 


Ibf 


B:= 


Mbf 


D:= 


D 


Mbfin 


Define V as the force and moment matrix. Force and moment matrix will be set to zero initially: 







Ny 


Ny 


0 

Fxy 


Nxy 


0 

Mx 


Mx 


0 

My 


My 


0 

.^Mxy j 


l,Mxy^ 


^0; 


Create a global matrix to represent the A, B, and D matrices: 

Global := stack(augment(A ,B) , augment(B ,D)) 



^ 266841.87 

49470.05 

47047.8 

0 

-0 

0 '' 


49470.05 

78650.66 

47047.8 

-0 

-0 

0 

Global = 

47047.8 

47047.8 

61638.11 

0 

0 

0 


0 

-0 

0 

12.29 

0.61 

0.39 


-0 

-0 

0 

0.61 

1.31 

0.39 


^ 0 

0 

0 

0.39 

o :39 

1.01 > 
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The global matrix is equivalent to: 


^Aii 

Ai 2 

Ai 3: 

Bn 

Bi 2 

Bi 3 

> 

to 

A-22 

^23; 

B21 

B22 

B23 

A31 

A32 

A33 

B31 

B32 

B33 

Bn 

Bi 2 

Bi 3 

Dll 

Di 2 

Di 3 

B21 

B22 

B23 

D2I 

D22 

D23 

vB 31 

B32 

B33 

D31 

D32 

D33 


Find the effective modulus in the x direction by assuming a unit force in the x direction: 


V = 




r Ny ^ 




X 


0 


Ny 


0 

Ibf 

Nxy 

:=V.M 

0 

in 

Mx 

in 

0 


My 




^Mxy J 



The Isolve function can be used to solve a linear system of n equations in n unknowns: 


^ SQx 


sOy 

TOxy 

Kx 

% 

V^xyy 


:= lsolve(Global,V) 


The Isolve function returns values for strain and curvature matrix. Young’s modulus and Poisson’s ratio 
can be calculated: 




Nx 

h-sox 


Ex 


= 1.13 X 10 


7M 
. 2 


in 
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Poisson’s ratio 




-spy 

5 px 


Vj^ = 3.17 X 10 


Now to find the effective modulus in the y direction, assume a unit force in the y direction: 


^ 0 ^ 
Ox 10 


V = 


1 X 10 


0 X 10 


Ox 10 


Ox 10 
VOX 10° y 


Ibf 


^ Nv ^ 




N 


xy 


Mv 


Mx; 


:= V- 


Ibf 


V^Mjjy J 


Now solve the system of equations for the strain matrix: 




spy 

TOxy 

Kx 

% 

V^xyy 


:= lsolve(Global,V) 


Young’s modulus in y direction: 




N„ 


ispy 


Ey = 2.0972 X 10 


6 Ibf 


Poisson’s ratio: 


2Qy _ 0 

Vyx:= Vy^:= 5.87x10 

spy 
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Now to find the effective modulus in the x direction, assume a unit force in the x-y direction: 


V = 





0 X 10 
0 X 10° 


% 

1 X 10° 

Ibf 

Nxy 

Ox 10° 

in 

Mx 

Ox 10° 


My 

,,0x 10°^ 


l^Mxy ^ 




Now solve the system of equations for the strain matrix: 


/'sn ^ 

sOx 


sOy 

70xy 

Kx 

K,. 


:= lsolve(Global,V) 


V^^xyy 


The shear modulus is calculated by: 


*^xy ■- 


Nxy 

70xy 


Gjto- = 3.22 X 10^ — 
^ . 2 
m 


D.2 Determination of Ply Strains 

Section D.l determined the effective material properties of the plate or laminate. It is also 
desirable to determine the stress and strain of each ply. This is accomplished by using the constitutive 
equations to determine the midplane strains and curvatures. 

For example, assume the laminated plate in the first example is loaded with 1,000 Ibf 
in the jc direction. The plate is 2 in. wide; the x direction is parallel to the 0° plies: 


Nv 


lOOOlbf 

Jin 


V, :« Nx — 
1 Ibf • 
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The force applied in the jc direction, set the first element of the V matrix equal to The values 
in the matrix must be unitless: 


/ 2^ 
2x 10 


Ox 10 


0 


Ox 10 


0 


0 X 10 


0 


Ox 10 


0 


„0 


VOx 10 y 

Now solve the system of equations for the strain matrix: 


^sox ^ 


sOy 

70xy 

Kx 

% 

V^^xyy 


:= lsolve(Global,V) 


Strains and curvatures at the midplane are found: 


EOx = 8 S2 X 10 

Kv = 0 X 10° 


-4 


sQy = -2.8 X 10 
Ky = 0 X 10° 


tOxy = -4-6 X 10 

Kxy=0x 10° 


-4 


Transform the midplane strains and curvatures into the strains in each ply. Now determine 
the strains in the first 45° ply. Zpjy is the distance from the midplane to the center of the ply: 


n ;= 2 




2ply 


;,= 2.5X 10 


-3 




50x ^ 


^Kx^ 

5y 


soy 

+ ^ly 

Ky 

Jxy> 


JOxyy 


l^KxyJ 
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8.82 :x 10 

5y 

= : 

-2.8 X 10“ ^ 

^Yxy> 


^-4.6 X 10“ 


Transform these strains into the principal material directions for the 45° ply: 
• Transformation matrix for 45° ply: 



f - 1 

- 1 

0 ^ 


5x 10 

5x10 

1 X 10 

T = 
n 

1 

o 

X 

1 

o 

X 

-1 X 10° 


- 1 

- 1 

0 


1,-5 X 10 

5x 10 

0x10^ 


• Transformation equation: 




f s > 
=x 




f 

7.13 X 10 

S2 

:= T • 
n 

5y 


82 

= 

5.31 X 10“ ^ 







^-5.81 X 10““^; 


• Use Reuter’s matrix to convert to engineering strain: 








f 

7.13 X 10 

S2 

:= R- 

B2 


B2 

= 

5.31 X 10“ ^ 

J12y 


^b12^ 


^712^ 


^-1.16 X 10“^; 


• Calculate the stresses within each ply. The stress can be calculated from the strain: 


i 

6 


6 


6 1 


6.56 X 

10 

4.55 X 

10 

4.7 X 

10 



6 


6 


6 

Ibf 

4.55 X 

10 

6.56 X 

10 

4.7 X 

10 

2 


6 


6 


6 

in 

^ 4.7 X 

10 

4.7 X 

10 

5.16 X 

10° J 
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Transformed reduced stiffness matrix for layer n = 2: 


^X 


f -4^ 

8.82 X 10 

By 



Jxy^ 


<-4.6k 10~^^ 


<^x 


Ex 

ay 

:= Qbar - 
^ n 

By 

VrxyJ 




• The stress for 45° ply in layer 2: 




'^x 




1 ^ 2.35 X 10^1 




1.99 X 10^ 

< 4.6 X 10^ ^ 


M 

. 2 
m 


• Transform the stress into the principal material direction: 




^ax'' 


<ai> 


f 3\ 

1.64 X 10 

a2 

:= T ■ 
n 

ay 


a2 

= 

7.23 X 10^ 

.^12 y 


i.^'i^xy j 




^,-1.16x 10^; 


D.3 Determine the Displacement of the Laminated Plate 

The definition of the plate curvatures provides an equation to determine the out-of-plane 
displacements for the plate. The curvature in the x and y directions is given by the following equations 


K3,:=-" w 
dx 


Ky:= 


dy 
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Since the curvatures are known values; i.e. constants, then w can be found by integration: 






■^2 


-Kvdx 


V / 


dx factor — ^ 


Example of the integration: 

• Limits of integration are along the edge of the plate: 


XI := 0 X2 := 0,.5.. 10 

yi ;= 0 72 := 0,.5..5 . 

• Example values ONLY! 

Kx := .1371* Ky := -.4733* Kjjy := 1 5* . 

• Integrate the a: curvature along the j edge of the plate: 


Wx(xi ,X2) := 


■X2 

r^2 




-Kxdx 


dx 


'Xl 


J 





^2 

Abiig Edge of Plate 
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Integrate the j curvature along the y edge of the plate: 


f 

■72. 

.-72 ■ 

'll 

\ 

-Kydy 





72 

Abng Edge of Plate 


The X and y deflections can be found from the strains and the curvatures: 


d d"' 

ex-— uo -2 — 
dx dx 


d d^ 

ey := — VQ -2 -w 

<ly dy 


_ d d 

Ixy := *2z — w— w 

dx dy 
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Midplane Strains 

- 4 

X 10 

- 4 

EQy = -2.S X 10 

- 4 

yOx7= JO 


Midplane Curvatures 

Kx = 0 X 10 ° 

Ky = 0 X 10 ° 

Kxy = 0 X 10° • 


• Remove the units: 

z 

^ l in 

Now integrating the midplane strains and curvatures over the length a: and the width j 
of the plate: 


Axfxi ,X2) 



^*2 



SQxdx-z^ 



-Kxdx 



w 

V 

xi ^ 




dx 


A)<0.0,10.Cr) =8.82x 10 


Deformed Shape of a Plate 



Abng Edge of Plate 
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Now integrating the midplane strains and curvatures for the j direction: 



"72 

r¥2 , 

rV2 

^y(yi .y2) - 

soy<^y-‘^^• 


-Kydy 

J 

yi 


J 


^yi 


Ay(0.0,5.IJ) =-1.4>{ 10 


Deformed Shape of a Plate 


a, 
• B 
Q 



Abng Edge of Plate 


D.4 Determination of the Nentral Axis 

The neutral axis is the plane of zero strain for any direction in an unsymmetric laminate: 
• Strain equations in matrix form: 




^ sQx ^ 



5y 

:= 

sOy 

+ z 

% 

Jxy> 


^70xy,j 
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Pure bending in the jc direction: 


f ° 1 


^Aii Ai2 Ai 3 Bii Bi2 B^^ 



0 


A21 A22 A 23 B21 B22 B23 


SOY 

0 


A31 A32 A33 B31 B32 B33 I 


tOxy 

Mx 


Bn Bi2 Bi 3 Dll D12 D13 


Kx 

0 


B21 B22 B23 D21 D22 D23 


% 

. 0 y 


^B3i B32 B33 D31 D32 D33J 


K^wJ 


• Solving for the midplane strains and curvatures: 


SQx '' 


''All Ai2 Ai 3 Bn Bi2 B13'' 

- 1 

f ° ^ 

SOY 


A21 A22 A23 B21 B22 B23 


0 

70xy 


A31 A32 A33 B31 B32 B33 


0 

Kx 


Bn Bi2 Bi 3 Dll D 12 D13 


Mx 

% 


B21 B22 B23 D21 D22 D23 


0 



^B 31 B32 B33 D31 D32 D33 ^ 


^ 0 2 


The equation can be reduced to the following: 


SQx ^ 


^BlP 

SOY 


B2I 

yOxy 


B3I 

Kx 


Dll 

% 


D2I 

Kxy J 


k° 3 l 2 


Invert the B and D matrices; if B is singular, the program will set = 0: 


Dmvrt D 


Binvit ■= B ™ ®rror B 



The above equation can then be separated to form an expression for the midplane strains and curvatures 


sOx 



■ 

<Kx> 



SOY 

= 


■Mx 

Ky 

= 

Dmvu^l 





l^KxyJ 




These equations are then substituted into the strain equation: 




^mvrt, , ' 




-]■ 

Sx 


1,1 



1,1 


Sy 

:= 

Bmvrt^ 1 

Mx + z 


Dmvrt^ 1 

■Mx 

Jxyy 








The general form of the equation to determine the neutral axis is given as guess values: 


^neuX •= 1 0 
^neuY ■= 1 Sy 

^neviX Y ■= 1 0 7xy 


:=0.0 

Mx:= 1.0 

:=0.0 

My := 0.0 

:=0.0 



The general form of the equations to determine the neutral axis is given as: 


^ X “ B invrt I I X B invrt I ^ ^ Y ^ ® invrt| ^ ^^xy - 

+ Z neuX (b invrt| |Mx+Dinvrt| ^My + Djnvrtj ^ ^ xy^ 

£ y “ B invrtj ^ M x + Binvrt^ ^ '^Y ^ ® '"'^2 3 ^ 

+ Z neuY ^ invrtj j x ^ Djnvrt^ 2 Y ^ ^ 3 ^ 

y xy = Binvrtj ^ M x + Bjnvrt^ ^ ^^y + B i„vrt^ ^ M xy ••• 

+ Z neuXT ^ M x + D jnvrt^ y + B invrt^ j '^xy^ • 

The neutral axis is found using the Find command: 


^ Z neuX 
Z neuY 


:— Find(Z neuX> Z neuY> Z neuXY) 


[ Z neuXY 
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The neutral axis for this example is at zero: 


^ neuX 


f 0 ^ 

2 neuY 

= 

-0 

^ Z neuXY ^ 


.0. 


Compare these results to the results from other methods: 
• Nettles equations: 


^ ^ 1,1 . "^ 1,2 
EKNettles := — ; — + — — 
h h 


1,3 


^2,3^1 ,3 "^1,2^3,3 

[ ^2 2 ^3 3- (A2 3) 


■^1,3 ^2,3'^1,2'^3,3 (^2,3) '^1,3 


3,3 A2^2'('^3,3) (^ 2 , 3 ) ■'^3,3 


ExNettles= 113x 

in 


^yWettles ■“ 


A A 

^^ 2,2 , ^^ 1,2 

’a a - a a 
^2,3^1,3 ^1,2 ^3,3 

h ' h 

r , 1 2 


^l,r^3,3" Pl,3) 


^,3 


■"^ 2,3 '^ 1 , 3 '^ 1 , 2 '^ 3,3 (^^1,3) '^ 2,3 


[ ^ l , r ['^ 3 , 3 ) ["^ 1 , 3 ) ■'^ 3,3 J 


EyNettles = 21 x 10 


6 Ibf 


Cl^tyNettlies 2 ' 


^3,3 (^2,3) ^'^1,3'^1,2'^2,2'^2,3 (^1,2) ['^2,3) [^1,3) '[^2,2) 


h h-A. 


2,2 


[^1,1) '[^ 2 , 2) [^1,2) '^ 2 , 2 _ 


GxyNettles ? 3.22 X 10^ — 
in 
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/ A . A \ 

^ , “1,3 “2,3 


vJ^Nettles := 




3,3 ; 


2,2 


(^2,3) 


3,3 


vxyNettles = -^ 34 x 10 


vyKNettks := 


*^^ 1 , 3 '^ 2,3 ^ 

A -2 

k ^ 3,3 J 


— -^1 
A 1 >1 

*3,3 


Halpin equations: 


vyxNettles = -5-27 x 10 


ExHalpin := 


'^ 1 , 1 '^ 2,2 ('^1,2) 
"^2,2'^ 


ExHalpm= 113 x 10 


7 M 

. 2 
in 


EyHalpin - 


^l,f^ 2,2 (^1,2) 


EyHalpin = 3^: K 10 


, 6 M 

.2 

m 


*33tyHalpin 
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184 




vxyHalpin - 

^ 2,2 

vxyHalpin = 62 ^ % 10 
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